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1.0  INTRODUCTION 


The  technique  of  sustaining  a  diffuse  discharge  by  electron  beam  Irradl- 

(1-3) 

atlon  of  a  gas  has  been  known  for  more  than  10  years.  In  this  type  of 

discharge,  a  high  energy  (>  100  kV)  electron  beam  Introduced  through  a  foil 
into  a  gas  maintains  a  controlled  level  of  Ionization  and,  hence,  a  controlled 
conductivity.  A  voltage  applied  between  electrodes  draws  a  current  through 
the  gas,  even  though  the  electric  field  Is  below  the  gas  breakdown  strength. 

The  e-beam-sustalned  discharge  has  a  number  of  advantages.  Because  the 
e-beam  current  and  the  discharge  voltage  can  be  Independently  controlled,  the 
technique  enables  one  to  Independently  adjust  the  electron  density  and  the 
electron  temperature  over  a  considerable  range.  Because,  unlike  a  conven¬ 
tional  glow  discharge,  this  discharge  Is  not  dependent  upon  wall  recombination 
for  stability.  It  can  be  scaled  to  large  volumes.  Unlike  a  preionized,  over- 
voltaged  avalanche  discharge,  the  e-beam-sustalned  system  Is  not  subject  to 
rapid  (<  1  ws)  arc  formation,  so  the  pulse  can  be  extended  for  tens,  or  even 
hundreds,  of  microseconds.  And,  since  a  high  energy  electron  beam  can  pene¬ 
trate  a  considerable  distance  through  gas  at  atmospheric  pressure  or  above, 
the  technique  can  be  used  to  produce  uniform,  diffuse  discharges  In  high 
pressure  gas.  There  Is  an  extensive  literature  on  e-beam-sustalned  dis¬ 
charges,  which  Is  sunmarlzed  in  several  reviews. 

The  e-beam-sustalned  discharge  has  been  used,  or  at  least  studied,  for 
many  applications,  including  gas  discharge  lasers,  discharge  chemical  process¬ 
ing,  the  measurement  of  gas  kinetic  processes,  and  as  an  electrical  switch 
that  can  be  opened  under  load.  Probably  the  most  extensively  explored  appli¬ 
cation  Is  the  CO^  Infrared  laser.  The  present  study  Is  based  on  a  prior 
Investigation  of  arcing  In  laser-type  discharges. But  the  use  of  an 
e-beam-sustalned  discharge  as  a  switch,  which  was  proposed  quite  early  In  the 
history  of  this  technology,  has  gained  Increasing  attention  In  recent  years. 

A  discussion  of  this  application,  with  references  to  the  earlier  work.  Is 

(9) 

given  In  the  recent  review  by  Kline.  The  purpose  of  the  present  study 
was  to  extend  the  Information  about  discharge  arcing  processes  reported  1n 
Refs.  7  and  B  to  the  somewhat  different  operating  regime  required  for  a  switch. 
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Switch  development  Is  of  Interest  because  a  compact,  efficient  opening 
switch  would  allow  one  to  store  energy  Inductively,  and  the  energy  density— 
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Joules  per  m  or  Joules  per  kilogram — that  can  be  stored  In  a  coll  is  more 
than  an  order  of  magnitude  greater  than  the  best  that  can  be  done  with  a 
capacitor  bank.^^^'^^^  An  opening  switch  would  also  be  useful  to  the  elec¬ 
tric  power  Industry  because  It  would  allow  the  Interruption  of  high  voltage  dc 
transmission  lines. 

An  e-beam-sustalned  discharge  can  be  used  as  an  opening  switch.  If  It 
Is  run  well  below  the  breakdown  voltage,  termination  of  the  e-beam  Interrupts 
the  discharge  current,  even  If  the  voltage  rises  because  of  the  inductive 
nature  of  the  storage.  Indeed,  for  pulsed  power  applications,  a  voltage  rise 
Is  desirable  since  the  stored  energy  can  then  be  delivered  to  the  load  at 
higher  voltage.  An  e-beam-sustalned  discharge  Is  always  somewhat  lossy  while 
conducting,  so  most  proposed  applications  envision  this  element  used  In  paral¬ 
lel  with  another  switch,  such  as  a  mechanical  circuit  breaker,  with  the  e-beam 
discharge  used  to  hold  down  the  switch  voltage  during  the  time  that  the  pri¬ 
mary  switch  Is  being  opened.  Unless  superconducting  colls  are  used.  Inductive 
storage  Is  also  dissipative,  and  would  be  used  only  as  a  stage  In  the  pulse¬ 
forming  system.  Typically,  one  might,  over  a  period  of  minutes,  accumulate 
energy  In  a  flywheel;  then.  In  a  time  -  1  sec,  transfer  this  energy  to  an  In¬ 
ductor;  and.  finally.  In  a  time  ~  1  tisec,  open-circuit  the  switch,  or  system 
of  switches,  forcing  the  Inductor  current  to  flow  through  the  load  (or  Into  a 
further  pulse  steepening  network  if  submlcrosecond  pulses  are  required). 

In  an  electric  power  dc  transmission  system,  the  e-beam-sustalned  switch 
would  also  be  used  In  parallel  with  a  circuit  breaker  to  conduct  the  current 
until  the  main  breaker  was  fully  open. 

The  requirements  for  a  switch  discharge  are  quite  different  from  a  laser 
discharge  since.  In  a  switch,  one  wants  to  minimize  the  energy  deposited  In 
the  discharge  medium.  The  gas  used  In  the  switch  should  have  a  high  breakdown 
strength  so  that  It  can  stand  off  the  voltage  after  the  switch  is  opened,  but 
also  have  a  high  conductivity  when  It  Is  e-beam-ionIzed  so  that  the  dissipa¬ 
tion  will  be  low.  One  way  to  achieve  the  latter  requirement  Is  to  use  a  gas 
with  a  Ramsauer  minimum  and  operate  the  switch  at  an  E/n  such  that  electron 
energies  are  1n  the  minimal  region  while  the  switch  Is  conducting.  To  further 
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reduce  the  collision  rate,  the  electron  energies  (l.e.,  the  electron  ran(:offl 
velocities)  should  be  low,  which  requires  a  gas  In  which  a  large  fraction  of 
the  electron  collisions  are  Inelastic.  For  this  cooling  process,  one  wants  a 
molecular  gas  In  which  there  are  vibrational  and  rotational  modes  that  can  be 
excited  by  low  energy  electrons.  Hence,  one  expects  that  a  molecular  gas  with 
a  Ramsauer  minimum  can  give  a  high  electron  drift  velocity  and,  hence,  a  high 
conductivity.  This  has  been  found  to  be  true  of  methane. and  most 
presently  envisioned  switch  designs  are  based  on  methane,  or  silane,  or  mix¬ 
tures  of  these  and  noble  gases. One  problem  with  such  gases  Is  that  a 
Ramsauer  minimum  can  make  the  dependence  of  electron  drift  velocity  upon  elec¬ 
tric  field  non-monotonIc,  which  can  make  the  discharge  unstable. The 
discharge  parameters  must  be  chosen  to  avoid  this  Instability. 

An  efficient  switch  would  be  recombination-dominated  during  the  current 
conduction  phase  but,  since  recombination  alone  gives  a  long  tall  to  the  Ioni¬ 
zation  after  e-beam  turnoff.  It  Is  usually  proposed  to  add  a  small  concentra¬ 
tion  of  an  attaching  species  to  terminate  more  efficiently  the  switch  conduc¬ 
tivity.  Since  the  electric  field  will  Increase  as  the  switch  Is  opened, It  has 
also  been  suggested  that  one  use  an  attacher,  attachment  rate  of  which  In¬ 
creases  with  C/n;  l.e.,  with  electron  energy.  (This  Is  true  of  dissociative 
attachment,  which  has  a  threshold  energy.)  Some  care  should  be  taken,  how¬ 
ever:  an  attachment  rate  that  Increases  with  C/n  can  make  the  discharge  un- 
(18) 

stable.  This,  however.  Is  less  of  a  concern  than  It  would  be  In  a  dis¬ 
charge  that  was  attachment-dominated  during  conduction. 

To  be  efficient,  the  switch  must  have  a  high  current  gain.  That  Is,  the 
discharge  current  must  be  greater  than  the  e-beam  current  needed  to  maintain 
the  conductivity.  In  a  recombination-dominated  discharge  at  a  given  gas  den¬ 
sity,  the  electron  density  and,  hence,  the  conductivity.  Increases  as  the 
square  root  of  the  e-beam  current.  This  favors  a  design  where  both  e-beam  and 
discharge  current  densities  are  low.  But  a  low  current  density  means  a  large 
switch  and.  If  the  switch  becomes  too  large,  the  advantage  of  Inductive  stor¬ 
age  Is  nullified.  An  absolute  limit  on  e-beam  current  and  pulse  length  Is  Im¬ 
posed  by  the  foil.  If  the  current  Is  too  high,  the  absorbed  energy  will  melt 
the  foil . 
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A  final,  absolute  requirement  is  that  the  switch  discharge  not  arc 
either  during  the  conduction  phase  or  afterwards,  when  the  voltage  increases 
as  the  switch  is  opened.  For  this,  it  is  obviously  necessary  that  the  peak 
voltage  be  kept  below  the  breakdown  threshold  of  the  gas.  But  that  is  not 
sufficient.  Although  it  is  run  below  the  breakdown  voltage,  an  e-beam  sus¬ 
tained  discharge  can  still  arc,  through  the  mechanism  of  discharge  streamers. 
These  streamers  are  hot,  filamentary  structures  that  start  at  the  electrodes 
and  propagate  across  the  discharge.  This  Is  the  process  that  was  studied  in 
our  earlier  work  on  laser-type  discharges.  The  same  phenomenon  has  also  been 
reported  by  others. 

In  discussions  of  e-beam  switch  capabilities,  two  alternative  simple 
criteria  are  sometimes  used  to  account  for  the  danger  of  arcing.  One  approach 
is  to  work  well  below  the  arcing  threshold;  i.e,  to  allow  a  large  "safety  fac¬ 
tor."  Unfortunately,  the  data  on  these  discharges  show  that  streamers  can 
still  form  and  that,  if  the  discharge  lasts  long  enough  and  deposits  suffi¬ 
cient  energy  in  the  gas,  it  can  still  arc. 

An  alternate  approach  is  to  place  a  limit  on  the  gas  heating,  noting 
that  laser-type  discharges  typically  arc  if  the  gas  is  heated  by  more  than 
~  500  K.  As  a  rough  estimate,  this  is  reasonable  but,  as  a  precise  design 
criterion,  it  is  inadequate.  As  discussed  in  Ref.  8,  the  specific  energy  that 
can  be  deposited  in  an  e-beam  sustained  discharge  before  arc  formation  depends 
upon  the  discharge  operating  conditions.  At  high  conductivity  (i.e.,  in 
higher  current,  lower  voltage  discharges),  the  allowable  gas  heating  is 
greater  than  for  laser  operating  conditions.  In  fact,  at  low  voltages,  the 
discharge  arcing  limit  was  found  to  be  Insensitive  to  discharge  power.  Such 
effects  have  also  been  reported  by  others. 

For  a  switch,  the  implications  are  good.  Not  only  does  the  higher  con¬ 
ductivity  operation  imply  that  less  energy  is  deposited  than  in  a  laser  dis¬ 
charge  with  the  same  current  density,  but  the  allowable  energy  deposition  may 
also  be  greater.  Note,  however,  that  this  advantage  is  somewhat  in  conflict 
with  the  need  for  high  current  gain.  The  latter  implies  that  low  e-beam  cur¬ 
rent  is  desirable,  while  the  need  for  a  compact,  arc-free  device  inclines  one 
towards  a  higher  e-beam  current.  It  is  clear  that  a  proper  assessment  must 
include  all  these  requirements  and,  in  particular,  must  be  based  on  a  more 
precise  understanding  of  the  arcing  limit. 
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2.0  SWITCH  SCALING 

Before  examining  the  arcing  mechanisms  in  more  detail,  however,  it  is 
useful  to  consider  the  switch  design  problem  more  systematically.  Suppose 
that  one  wishes  to  store  a  megajoule  and  then  deliver  this  energy  to  a  load  at 
200  kV  in  10  ys.  If  a  simple  Inductor  is  the  store  and  the  load  voltage  is 
constant,  the  current  would  ramp  down,  starting  at  one  megampere  Just  after 
the  switch  opens.  So  the  switch  must  conduct  lO^A,  with  the  voltage,  of 
course,  much  less  than  200  kV  while  the  switch  is  dosed.  Note  that  the  vol- 
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ume  of  the  storage  coil  can  be  as  small  as  0.1  m  ;  a  capacitor  bank  that  can 

6  3 

store  10  J  would  occupy  at  least  1  m  . 

The  e-beam-sustained  switch  must  conduct  for  a  time  sufficient  to  open  a 
parallel  breaker,  then  open  in  a  time  much  shorter  than  the  time  to  deliver 
the  energy  to  the  load.  Since  we  are  postulating  a  10-tis  load  delivery,  we 
assume  that  the  switch  must  open  in  1  \ts.  To  open  a  parallel  breaker,  one 
would  like  a  conduction  time  of  a  few  ms  but,  since  this  is  probably  imprac¬ 
tical,  we  require  that  the  e-beam  switch  conduct  for  100  ys,  which  is  at  least 
long  enough  to  open  an  explosive  switch, and  perhaps  to  blow  out  an  arc 
or  open  a  fast  mechanical  device. 

Assume  that  the  switch  is  to  run  on  90  percent  CH^,  10  percent  A  and, 
for  efficiency,  is  to  be  recombination-dominated  during  the  conduction  phase. 
Then,  just  after  the  e-beam  ceases,  the  electron  density  in  the  discharge 
evolves  according  to 


dn 


n„  dt 
e 


-an. 


(1) 


where  a  is  the  recombination  rate  coefficient.  Since  the  switch  must  open  in 
1  ys,  we  require  that  the  initial  rate  of  decay  be  at  least  2  x  10^/sec. 

This  limits  the  discharge  electron  density;  >  2  x  10^ /a. 

The  recombination  coefficient,  a,  for  methane  has  been  calculated  to  be 

-7  3  (9)  (20) 

~  2  X  10  cm  /sec,  but  recent  measurements  by  Bletzinger  gave 
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a  value  an  order  of  magnitude  larger.  Using  this  measured  value, 

-6  3 

2  X  10  cm  /sec,  one  concludes  that  the  electron  density  must  be 

12  3 

n  >  10  electrons/cm  . 
e  ~ 

Note  that  this  determines  the  e-beam  power  deposition  in  the  discharge. 

12 

If  the  electron  density  is  10  and  the  plasma  is  recombining  at  a  rate 

6  18  3 

d(ln  n^l/dt  =  -  2  x  10  /sec,  one  must  produce  2  x  10  electrons/cm  sec 

to  keep  the  switch  closed.  If  E^,  the  effective  ionization  potential  (the 

e-beam  energy  deposited  per  electron  produced)  is  35  eV  (a  typical  number), 

3 

then  the  e-beam  must  deposit  11.2  watts/cra  while  the  switch  is  closed. 

If  the  switch  is  to  run  at  an  E/n  which  puts  electrons  in  the  Ramsauer 
minimum  of  the  gas.  the  discharge  current,  j^,  is  also  fixed  by  n^.  For 
methane-argon,  9:1,  the  minimum  starts  at  an  E/n  of  2.5  Td  and  gives  an  elec¬ 
tron  drift  velocity  of  approximately  10^  cm/sec. Hence,  the  discharge 

2 

current  density  is  j-  =  e  n  v  =1.6  A/cm  .  This,  in  turn,  defines 

^  ^  ^  6  2 
the  area  of  the  switch  needed  to  carry  10  A:  the  discharge  must  cover  62.5  m  . 

This  is  a  major  disadvantage.  A  switch  this  size  would  completely  null¬ 
ify  the  advantages  of  compact  storage  offered  by  inductors. 

What  can  be  done  to  improve  the  situation?  There  are  two  changes  that 
will  help.  The  first  change  is  to  deliver  the  energy  at  higher  voltage  and 
lower  current.  The  second  change  one  can  effect  is  to  run  the  discharge  at  a 
higher  electron  density  to  increase  the  current.  In  practice,  one  would 

almost  certainly  have  to  do  both  to  make  the  switch  small  enough  to  be  useful. 

6 

It  is  clear  that  the  circuit  parameters  assumed  above  (10  A  while 

closed,  200  kV  after  opening)  are  not  appropriate  to  this  type  of  energy  stor- 

5 

age.  Suppose  instead  that  we  specify  2  x  10  A  conduction  current  and  1  MV 

open  circuit  standoff.  Suppose  also  that  the  switch  is  run  at  an  electron 
12  3 

density  of  5  x  10  /cm  .  Then  the  e-beam  power  deposition  (assuming 
-6  3 

a  =  2  X  10  )  is  280  W/cm  ,  the  discharge  current  density  (at  E/n  »  2.5  Td) 

2  2 
is  8  A/cm  ,  and  the  area  of  the  switch  is  2.5  m  .  Since  one  reason  for 

3 

using  this  switch  is  to  replace  a  few  m  of  capacitors  by  a  smaller  induc¬ 
tor,  this  switch  is  still  undesirably  large,  but  at  least  we  have  reduced  it 
to  a  size  that  begins  to  be  practical. 

The  electric  field,  the  discharge  electrode  spacing,  the  discharge  power 
density,  and  the  required  e-beam  current  density  all  depend  upon  the  switch 
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gas  density.  At  the  desired  E/n,  the  electric  field  is  625  p/p  V/cm,  where 

0 

p/p^  is  the  ratio  of  switch  gas  density  to  the  density  at  atmospheric  condi¬ 
tions.  The  discharge  power  dissipation  is 


P  =  Ej  =  5  p/p^  kW/cm^  (2) 

For  a  high  pressure  discharge,  this  is  much  larger  than  the  e-beam  power 
deposition.  Using  the  simple  arcing  criterion  mentioned  earlier  -  that  the 
gas  should  not  be  heated  by  more  than  500°K  -  gives,  for  typical  molecular 
gases,  an  allowable  energy  deposition, 

W  <  0.6  p/p^  J/cm^  (3) 

A  significant  amount  of  energy  could  be  deposited  when  the  switch  opens, 

if  the  voltage  rises  while  the  switch  is  still  conducting.  Indeed,  it  can  be 
(22) 

shown  that  the  external  circuit  imposes  an  irreducible  switch  dissipa¬ 
tion  during  opening.  But,  because  the  heating  of  the  gas  is  a  major  design 
constraint  and  because  heating  during  opening  would  be  at  high  E,  low  J,  where 
the  allowable  heating  is  less,  it  is  unlikely  that  a  large  dissipation  upon 
opening  would  be  acceptable.  It  is  more  likely  that  the  external  circuit 
would  be  designed  to  prevent  such  heating.  One  could,  for  example,  put  an 
R-C,  or  a  more  complex  waveform-shaping  network,  in  parallel  with  the  switch 
to  hold  down  the  voltage  for  ~  1  ns  while  the  switch  was  opening.  (If  this 
unacceptably  degraded  the  load  voltage  rise  time,  a  spark  gap  or  other  closing 
switch  could  be  added  between  this  system  and  the  load.)  This  requires  some 
capacitive  energy  storage  in  the  filter,  but  the  stored  energy  would  be  only 
~  0.1  of  the  total  in  the  coil,  so  the  filter  capacitor  would  still  be  much 
smaller  than  the  capacitor  that  would  be  needed  to  store  all  the  energy. 

Hence  we  assume  that  the  significant  dissipation  in  the  switch  occurs  during 
the  conduction  phase. 

The  allowable  energy  dissipation  and  the  switch  power  together  put  a 
limit  on  the  time  the  switch  can  remain  closed  without  arcing.  For  our 
example, 


t 


W  , 

p  <  120  ns 


(4) 
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Note  that  this  time  Is  Independent  of  the  gas  density.  This  pulse  length 

meets  the  requirement  of  100  pS  assumed  earlier,  but  not  by  much.  One  would 

like  a  longer  pulse  and,  If  the  arcing  limit  can  be  Improved,  that  could  be 

done.  The  overall  scaling  should  be  noted:  since  the  electric  field  and  the 

allowable  energy  deposition  are  both  proportional  to  gas  density,  while  the 

discharge  current  density,  at  fixed  E/n,  Is  proportional  to  n  ,  the  allowed 

e 

pulse  length  varies  Inversely  with  the  discharge  electron  density.  This  Is 
the  essential  tradeoff:  by  Increasing  n^  at  a  given  gas  density  one  obtains 
a  more  compa'.:t  device  (and  a  more  rapid  opening),  but  at  the  price  of  a  reduc¬ 
tion  in  tho  time  the  switch  can  remain  closed. 

The  e-beam  current  needed  to  maintain  n  Is  found  from  the  equation 

e 


a  n^  = 


*eb 


e  E, 


dm 


(5) 


6  2 

where  dE/dm  Is  the  mass  stopping  power  of  the  gas  -  typically  2  x  10  eV  cm  /g, 
p  Is  gas  density,  and  E.  Is  the  effective  Ionization  potential  used 
earlier.  Taking  E^  =  35  eV  and  a  =  2x10  cm  /sec,  one  has 

2 

J,|,  “  6.7x10-2’ 

12 

For  Og  »  5x10  ,  this  gives 


The  current  gain  Is 


167 


'eb 


(P/Po) 


mA/cm 


(7) 


e  V  n 
e  e 


^eb  6.7  X  lO’^^  n2/(p/p^) 


2.4  X  10 
n_ 
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(8) 


For  the  assumed  n 


e’ 


7^  =  48  (p/p  ) 
Jeb 


(9) 


At  10  atmospheres,  the  current  gain  Is  almost  500.  To  Improve  the  gain,  one 

can  either  reduce  n  -  I.e.,  work  at  low  j  .  and  low  j.  -,  or  Increase 

e  ep  D 

the  density.  Since  lowering  n  makes  the  switch  Impractical ly  large,  one 

e 

wants  to  work  at  higher  density.  The  same  Is  true  of  the  foil-melting  limit. 
12 

Ifn^-SxlO  Is  needed  to  make  a  compact  system  or  to  reduce  the 

switch-opening  time,  the  foil-melting  limit  (which  gives  a  maximum  J  .  for  a 

60 

given  pulse  duration)  then  Implies  a  minimum  p/p  for  the  system.  High 

0 

pressure  operation,  of  course,  requires  a  strong  foil  to  hold  the  gas  without 

(23) 

breaking,  but  this  can  be  done.  A  recent  study  showed  that  more  than  35 
atmospheres  can  be  held  by  an  0.6  mil  T1  alloy  foil  on  a  frame  with  0.3  cm 
slots.  ' 


If'Jeb  is  kept  constant  (e.g.,  held  at  the  foil  melting  limit),  the 


electron  density  scales  ast/’pTp”,  so  the  recombination  (switch-opening) 
time  is  Inversely  proportional  to^  p/p^  and.  If  E/n  Is  kept  constant  (I.e., 
If  the  gap  spacing  Is  decreased  Inversely  with  p),  the  current  density, 
varies  as  y^p/p^,  so  the  current  gain  Increases  as  V^p/p^  and  the  switch 
area  varies  Inversely  with  this  square  root  of  the  density.  It  Is  clear  that 
the  switch  should  be  designed  to  use  the  highest  practical  pressure. 

The  gap  spacing  between  switch  electrodes,  which,  for  a  given  p  and  E/n 
determines  both  the  conduction  voltage  of  the  switch  and  the  e-beam  voltage 
needed  to  maintain  a  nearly  uniform  Ionization,  Is  set  by  the  electric  field 
which  the  switch  can  hold  off  after  opening.  The  Paschen  curve  for  methane 

Assuming  that,  after  a 


4  (9) 

shows  a  breakdown  field  of  2.5  x  10  V/cm  atm. 

4 

pulse,  one  can  do  half  this  well,  1.25  x  10  gives  the  gap  spacing  for  a 


1-MV  system: 


So,  for  10  atmospheres,  the  anode-cathode  spacing  Is  8  cm.  To  penetrate  this 
gas,  the  e-beam  must  have  an  energy  of  -  250  keV.  This  Is  practical.  One 
does  not,  however,  want  to  use  a  much  higher  voltage.  An  e-beam  with  a  few 
times  250  kV  would  produce  hard  X-rays,  which  would  create  a  shielding  prob¬ 
lem.  Hence  the  discharge  voltage  assumed  here  (1  MV  open  circuit)  Is  an 
attractive  design  point. 

For  this  gap  spacing  and  E/n,  the  closed  circuit  voltage  Is  50  kV,  the 

3  12 

switch  volume,  A.d,  Is  2.0  (p./p)  m  (at  constant  n.  >  5  x  10  ),  the 

®  10  ® 

total  discharge  power,  I.V  >  1.0  x  10  W,  and,  for  a  100  vs  pulse,  the  dis¬ 
sipated  energy,  I.V.t  -  1.0  NJ,  equal  to  the  energy  to  be  delivered  to  the 
load. 

Note  that  this  last  result  Is  Independent  of  most  of  the  switch  para¬ 
meters.  For  a  given  ratio  of  open  E  to  closed  E  In  the  switch,  and  a  given 
ratio  of  switch  conduction  time,  t,  to  load  delivery  time,  t^^  (5  ps  In  our 
example,  since  the  current  ramps  down  In  10  t>s),  one  has  a  system  *Q* 


In  fact,  the  efficiency  Is  even  worse  than  this.  Since  one  Initially 
must  store  2  HJ,  the  coll  —  and  hence  switch  —  current  Is  Initially  greater 
by  2  and  the  Initial  dissipation  will  be  higher.  Note  that  this  also  reduces 
the  time  that  the  switch  can  conduct  without  arcing. 

For  applications.  It  may  be  acceptable  to  dissipate  half  the  energy  In 
the  switch.  If  the  cost  of  spinning  up  a  flywheel  and  then  storing  that  energy 
In  a  coll  Is  relatively  small,  but  a  higher  efficiency  Is  certainly  prefer¬ 
able.  The  above  simple  result  shows  that  the  efficiency  can  be  Improved  only 
by  reducing  the  ratio  of  conduction  time  to  load  delivery  time,  or  else  by  In¬ 
creasing  the  ratio  of  open-to-closed  circuit  voltage.  The  ratio  of  times  Is 
already  Inconvenient.  One  would  like  the  switch  to  conduct  for  a  longer  time, 
and  It  may  be  necessary  to  deliver  the  energy  to  the  load  In  a  shorter  time. 


10 


A  better  change  Is  to  Increase  E  /E  ,  This  may  be  possible.  When 

open  closed 

open,  the  switch  may  be  able  to  hold  off  more  than  half  the  breakdown  vol¬ 
tage.  More  Important,  during  conduction.  It  may  be  practical  to  use  a  lower  E 
than  assumed  here.  The  plots  of  electron  drift  versus  E/n  In  Ref.  20  have  a 
plateau  (actually  a  slight  decline)  above  2.5  Td,  so  this  value  was  taken  as 
the  "Ramsauer  minimum,"  the  point  where  a  switch  should  run.  But  below 
2.5  Td,  the  curve  looks  linear.  The  drift  velocity  Is  just  proportional 
to  E.  If  true,  this  means  that  the  gas  conductivity  Is  constant  there.  In  a 
sense,  2.5  Td  Is  the  upper  bound  of  the  Ramsauer  minimum.  At  higher  fields, 
the  drift  velocity  does  not  Increase  because  electrons  have  energies  higher 
than  the  minimum  and  hence  are  scattered  more  efficiently.  So  2.5  Td  Is  the 
highest  field  that  one  should  consider  for  the  conduction  phase.  To  use  a 
lower  field,  one  must  of  course  accept  a  lower  and,  hence,  a  larger 

switch,  but  the  tradeoff  Is  an  Increase  In  E  /E  ,  I.e.,  an  In¬ 

open  closed 

crease  In  the  switch  efficiency. 

The  simple  coll  storage  system  assumed  here  may  also  not  be  the  best 
choice,  since  the  triangular  load  delivery  waveform  Is  unlikely  to  be  what  one 
wants.  It  may  be  better  to  use  a  lumped  element  transmission  line  In  which 
the  energy  Initially  Is  stored  Inductively.  Another  possibility  Is  to  store 
more  energy  In  a  larger  coll.  Then  the  current  would  not  droop  significantly 
during  either  the  conduction  time  of  the  switch  or  the  load  delivery  Interval, 
and  one  might  be  able  to  return  the  remainder  to  the  flywheel  or  use  It  on 
successive  pulses  of  the  system.  Clearly,  switch  development  and  more 
detailed  work  on  the  circuits  In  which  such  switches  would  be  used  should  be 
done  together. 

Therefore,  the  foregoing  results  for  a  simple  coll  circuit  are  not 
final,  but  they  are  useful  as  a  guide  to  the  Improvements  needed  for  a  prac¬ 
tical  switch.  The  principal  limitations  Indicated  by  this  analysis  are  the 
physical  size  and  the  energy  efficiency  of  the  switch.  To  reduce  the  size  re¬ 
quires  a  high  current  density,  which  Is  best  achieved  by  a  high  gas  density  l\t 
the  device.  The  device  also  becomes  more  compact  as  one  goes  to  higher  vol¬ 
tage  and  lower  current.  But  the  useable  open  circuit  voltage  Is  limited  by 
the  need  to  keep  the  e-beam  voltage  below  the  level  at  which  X-ray  shielding 
problems  would  become  a  major  consideration. 
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The  other  major  limitation  Is  the  switch  conduction  time,  which  must  be 
long  enough  to  open  a  parallel,  lower-resistance  switch.  Both  the  open- 
circuit  voltage  and  the  conduction  time  are  limited  by  streamer  growth,  and 
there  Is  probably  a  tradeoff:  If  the  switch  conducts  for  a  shorter  time, 
streamers  will  not  grow  as  far,  and  the  switch  can  probably  then  hold  off  more 
voltage  after  opening.  The  choice  of  closed  circuit  voltage  (I.e.,  E/n) 
affects  all  three  parameters:  If  the  switch  Is  run  at  a  lower  E/n,  the  area 
must  be  larger,  but  the  conduction  time,  or  the  overall  system  efficiency,  or 
both,  can  be  Increased. 

Switch-opening  time,  on  the  other  hand,  does  not  appear  to  be  a  problem, 
at  least  If  1  gs  Is  adequate.  At  the  high  gas  density  and  relatively  high 
electron  density  needed  for  a  compact,  efficient  system,  the  recombination 
rate  Is  rapid.  A  small  amount  of  attacher  may  be  needed  to  truncate  the  re¬ 
combination  tall,  but  a  larger  quantity,  which  would  accelerate  the  Initial 
turnoff,  is  unnecessary. 

The  strongest  conclusion  of  these  design  considerations  Is  the  need  to 
go  to  higher  pressures.  Three  major  characteristics  that  a  practical  switch 
must  have  (high  current  density,  high  current  gain,  and  rapid  opening)  are  all 
Improved  by  Increasing  the  density  of  the  discharge  gas. 

The  essential  unanswered  questions  that  will  determine  the  practicality 
of  such  a  switch  are:  the  realistic  limit  on  the  pressure  at  which  the  switch 
can  be  run;  and  the  arcing  limit,  the  energy  deposition  and  the  open  circuit 
E/n  which  the  switch  can  sustain  without  arcing.  Of  particular  Importance  Is 
the  pressure-scaling  of  this  limit.  The  gains  apparently  offered  by  an  In¬ 
crease  In  the  pressure  are  real  only  If  the  specific  energy  deposition  can  be 
maintained  at  higher  densities. 
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3.0  THE  STREAMER  INSTABILITY 


3.1  EXPERIMENTAL  RESULTS  AND  DISCUSSION 

The  theoretical  description  of  arc  formation  In  an  e-beam-sustalned  dis¬ 
charge  developed  under  this  program  Is  based  on  observations  of  streamer  form¬ 
ation  In  a  small  discharge  device.  That  work  Is  described  In  Refs.  7  and  8 
but,  for  completeness.  It  will  also  be  summarized  here. 

High-speed  framing  camera  photographs  of  the  discharge  showed  that  arc¬ 
ing  results  from  discharge  streamers,  luminous  filaments  that  start  at  the 
electrodes  (usually  the  cathode,  although  anode  streamers  can  also  be  pro¬ 
duced)  and  grow  Into  the  discharge.  An  example  of  such  a  sequence  of  photos 
Is  shown  In  Figure  1.  The  same  phenomenon  had  been  seen  earlier  In  open  shut¬ 
ter  photographs  of  pulsed  discharges  when  the  pulse  length  was  short  enough 
that  the  discharge  didn't  arc.  Streamers  In  e-beam-sustalned  discharges  have 
also  been  observed  by  others. 

As  can  be  seen  In  Figure  1,  the  streamer  has  a  characteristic  struc¬ 
ture.  It  consists  of  a  filamentary  core  with  a  glowing  ball  on  the  end.  In 
color  photographs  (shown  In  Ref.  8),  the  core  Is  seen  to  be  brilliant  white, 
resembling  an  arc,  while  the  halo,  or  cap  sphere.  Is,  at  least  In  nitrogen,  a 
reddish  or  salmon  color,  more  closely  resembling  a  corona-type  of  discharge. 

As  the  streamer  grows,  the  form  remains  the  same,  while  the  column  and  the 
halo  are  seen  to  grow  proportionately. 

When  the  length  of  the  streamer  Is  plotted  against  time,  as  In  Figure  2. 
the  growth  Is  found  to  be  exponential.  This  was  always  the  case  In  our  exper¬ 
iments,  except  In  a  few  high  current  discharges  where  the  voltage  drooped  dur¬ 
ing  the  pulse.  Typically,  the  discharge  arcs  when  the  streamer  has  grown 
about  two  thirds  of  the  way  across  the  gap.  The  evident  explanation  is  that 
the  streamer  Is  a  highly  conducting  structure,  essentially  at  cathode  poten¬ 
tial,  and  an  arc  forms  when  the  field  across  the  remaining  gap  between  the 
streamer  and  the  anode  exceeds  the  threshold  for  rapid,  avalanche-type  break¬ 
downs. 
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Figure  1  Streamer  growth.  The  marks  on  the  right  denote  successive  frames. 
The  framing  rate  is  2?, 800  frames/sec.  The  line  of  sight  Is 
parallel  to  the  cathode  surface.  The  anode  to  cathode  spacing  was 
1  Inch  and  the  view  here  includes  most  of  the  discharge. 


LENGTH  (cm) 


So  streamer  growth  Is  characterized  by  a  growth  rate,  or  exponentiation 

time,  rather  than  by  a  velocity.  This  behavior  Is  different  from  the  arcing 

of  an  overvoltaged,  avalanche  discharge,  where  the  filamentary  streamers  that 

precede  an  arc  seem  to  advance  at  a  roughly  constant,  although  much  higher, 
(251 

velocity.  In  our  studies  of  the  slower  streamer  propagation  in  e-beam- 
sustained  discharges,  no  velocity  effects  were  seen,  even  for  streamers  that 
became  supersonic  with  respect  to  the  background  gas.  There  was  no  change  in 
the  exponential  growth  and.  in  shadowgraphs,  which  were  also  taken  with  the 
framing  camera,  there  was  no  evidence  of  shock  waves.  Apparently  the  gas,  at 
least  near  the  streamer,  was  heated  enough  so  that  the  gas  motion  was  still 
subsonic.  As  discussed  below,  streamers  in  the  experiments  were  seen  to  form 
rapidly — in  the  first  few  microseconds — with  a  size  of  about  1  millimeter,  and 
the  discharge  arcing  time  was  about  2.5  times  the  streamer  exponentiation  time 
seen  in  the  photographic  sequences. 

The  streamer  growth  rate  was  seen  to  depend  upon  both  the  voltage  and 
the  power  of  the  discharge,  as  shown  in  Figure  3,  where  data  in  three  differ¬ 
ent  growth  rate  ranges  are  plotted  against  power  and  electric  field.  The 
solid  line  is  a  contour  of  constant  growth  rate,  as  determined  from  the  data. 
What  this  shows  Is  that  one  cannot  avoid  arcing  simply  by  limiting  the  elec¬ 
tric  field  or  the  gas  heating.  The  correct  criterion  Involves  both.  At  low  E 
field,  one  can  allow  more  gas  heating.  Indeed,  at  low  E,  where  a  switch  would 
operate  while  closed,  the  line  in  Figure  3  is  nearly  vertical.  Implying  that 
in  this  regime  the  streamer  growth  rate  is  largely  Independent  of  discharge 
power.  (Experimentally,  at  constant  E  the  power  is  varied  by  changing  the 
e-beam  current,  which  changes  the  discharge  gas  Ionization,  and  hence  the  con¬ 
ductivity,  and  hence  the  discharge  current.)  At  high  E,  where  laser  dis¬ 
charges  typically  run,  the  line  is  nearly  horizontal,  implying  that  the 
streamer  growth  rate  Is  primarily  power-dependent  In  that  regime.  As  noted 
above,  such  variation  and.  In  particular,  the  fact  that  more  gas  heating  Is 
allowable  at  lower  E,  has  also  been  seen  by  other  Investigators. 

The  evident  explanation  is  that  the  streamer  propagation  Involves  both 
field  ionization — i.e.,  avalanching — and  gas  heating,  and  that  the  growth  rate 
is  most  strongly  dependent  upon  whichever  process  Is  weaker.  To  be  useful,  a 
theoretical  model  must  explain  this  behavior  and  predict  the  variation  with 
pressure,  gas  composition,  etc. 
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r  ,  the  streamer  growth  rate,  is  in  units 

OF  (100  8)-1  THE  ALLOWABLE  POWER  DEPOSI¬ 
TION  IS  HIGHER  AT  LOW  E/N.  THIS  IMPLIES  THAT 
THE  SWITCH  IS  MOST  LIKELY  TO  ARC  DURING 
THE  OPENING  PHASE. 


Figure  3.  Dependence  of  streamer  growth  rate  upon  discharge  power  density  and 
electric  field,  y  1s  In  units  of  (100  vs)-1. 
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Streamers,  once  formed,  can  be  maintained,  and  can  continue  to  grow 
after  the  current  Is  almost  completely  shut  off.  In  the  experiments,  post 
arcs  were  seen  to  be  due  to  streamers  that  continued  growing  slowly  after  the 
end  of  the  discharge.  (In  this  device,  the  main  discharge  ended  when  the 
e-beam  pulse  ended,  but  the  discharge  voltage  remained  on  for  a  considerable 
time  afterwards.)  An  example  Is  shown  In  Figure  4,  a  framing  sequence  of  a 
discharge  that  arced  about  a  millisecond  after  the  end  of  the  e-beam  pulse. 

At  these  voltages,  the  discharge  never  pre-arced,  even  though  the  voltage, 
applied  by  a  mechanical  switch,  was  present  for  several  seconds  before  the 
pulse.  The  post  arc  was  clearly  caused  by  the  discharge  and,  as  can  be  seen, 
was  due  to  a  streamer,  a  phenomenon  that  must  be  considered  If  an  e-beam  sus¬ 
tained  discharge  Is  to  be  used  as  a  repetitively  pulsed  switch. 

Qualitatively,  we  think  that  a  streamer  Is  a  composite  structure  around 
a  thermally  Ionized  core,  essentially  an  arc,  which  grows  Into  the  discharge 
by  heating  the  gas  ahead  of  It.  Because  the  streamer  column  is  a  good  conduc¬ 
tor,  It  collects  a  concentrated  current,  as  would  a  wire  stuck  Into  the  dis¬ 
charge  from  the  cathode.  Indeed,  In  some  of  the  early  experiments,  a  wire 
loop  about  1  mm  In  diameter  was  mounted  on  the  cathode  and  the  largest 
streamers  always  grew  from  there,  which  simplified  the  observations. 

The  concentrated  current  through  the  streamer  maintains  the  temperature 
and  hence  conductivity  of  the  core  and  also  heats  the  gas  around  the  tip. 
Because  the  streamer  Is  a  good  conductor.  It  Is  close  to  cathode  potential,  so 
there  Is  an  enhanced  electric  field  around  the  tip.  One  expects  this  field  to 
produce  non-thermal  Ionization,  which  Is  consistent  with  the  glowing  colored 
ball  seen  In  the  photographs.  In  this  region  one  expects  the  field  to  be 
close  to  the  gas  breakdown  strength.  A  higher  field  would  give  rapid  buildup 
of  the  Ionization,  which  would  act  to  reduce  the  field;  a  lower  field  would 
allow  this  non-thermal  Ionization  to  decay,  which  would  cause  the  field  to  In¬ 
crease.  Since  the  voltage  drop  around  the  streamer  tip  1s  proportional  to  the 
length  of  the  streamer  (It  Is  Just  the  voltage  that  would  have  been  dropped 
uniformly  between  that  point  and  the  cathode  If  there  were  no  streamer),  a 
clamped  electric  field  in  the  halo  Implies  that  the  halo  radius  should  be  pro¬ 
portional  to  streamer  length,  which  agrees  with  the  observations. 


Note  that  this  does  not  imply  a  constant  E  within  the  halo,  since  the 
gas  temperature  certainly  varies  there,  and  the  field  needed  for  avalanching 
is  proportional  to  density.  But  if  the  temperature  profile  is  self-similar  as 
the  streamer  grows  (and  we  think  it  is),  the  E  field  profile  will  be,  too. 

It  is  plausible  that  such  a  structure  would  grow  exponentially.  For 
most  of  its  growth  time  the  whole  streamer  is  relatively  close  to  the  cathode, 
so  it  is  effectively  growing  into  a  semi-infinite  discharge.  A  good  conduc¬ 
tor,  the  streamer  will  collect  roughly  that  current  that  would  otherwise  have 
gone  to  the  area  on  the  cathode  that  is  "shadowed"  by  the  streamer.  That  is, 
it  will  collect  a  current  roughly  equal  to  the  background  current  density 
times  an  area  with  radius  equal  to  the  length  of  the  streamer.  Hence  the 
streamer  current  grows  as  the  square  of  the  length  of  the  streamer. 

Since  the  radius  of  the  cap  sphere,  or  halo,  is  also  proportional  to 
streamer  length,  its  surface  area  also  grows  as  the  streamer  length  squared, 
so  the  current  density  there — and  at  other  regions  within  and  around  the 
halo — remains  constant  or,  more  precisely,  self-similar  as  the  streamer 
grows.  Since  this  is  also  true  of  the  electric  field,  it  is  also  true  of  the 
power,  or  gas  heating.  If  the  heating  profile  within  and  around  the  halo  is 
self similar  as  the  streamer  grows,  the  time  to  heat  the  gas  In  the  halo  or,  at 
least,  in  the  leading  part  of  the  halo  up  to  the  core  temperature.  Is  constant 
as  the  streamer  grows.  In  each  such  heating  time,  the  streamer  grows  by  a 
halo  radius,  a  distance  proportional  to  streamer  length.  A  structure  that 
grows  by  a  certain  fraction  of  its  length  in  a  fixed  time,  of  course,  grows 
exponentially,  as  streamers  are  consistently  observed  to  do.  A  proper  theo¬ 
retical  analysis  should  test  and  quantify  this  picture. 

3.2  STREAMER  INITIATION 

In  short-pulse,  high-voltage  discharges,  such  as  are  used  in  TEA  lasers, 
careful  shaping  and  polishing  of  the  electrodes  are  crucial  to  the  prevention 
of  premature  arcing  of  the  discharge.  This  suggests  that  similar  techniques 
might  be  effective  in  preventing  or  at  least  delaying  the  formation  of 
streamers  in  e-beam-sustained  discharges.  However,  a  more  careful  analysis  of 
the  problem,  plus  some  observations  made  in  our  experimental  work,  make  this 
solution  seem  less  promising. 
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The  cathodes  used  in  our  experiments  were  not  carefully  shaped  for  field 
uniformity,  and  the  field  was  highest  at  the  edges.  Yet,  as  shown  by  burn 
spots  left  by  arcs,  streamers  did  not  always  start  at  the  edges.  Many  of  them 
did  but,  as  can  be  seen  from  Figure  5,  a  photograph  of  a  cathode  after  a  ser¬ 
ies  of  discharges,  there  were  also  burn  marks  scattered  over  the  flat  central 
region  where  the  cathode  E  field  (or  at  least  the  vacuum  field)  would  have 
been  weakest.  The  cathode  shows  a  tendency  for  the  arcs  to  clump  in  the  cen¬ 
ter,  where  the  e-beam  and,  hence,  the  discharge  current  density,  was  highest. 

Consideration  of  the  structure  and  evolution  of  the  cathode  layer  in  a 
high-pressure  discharge  suggests  two  mechanisms  by  which  streamers  might  be 
started:  growth  from  cathode  spots  and  an  instability  of  a  large-area  cathode 
layer.  There  is  evidence  of  both  mechanisms  in  the  experimental  data.  Work 
on  low  pressure  discharges'  ’  ’  has  provided  a  fairly  complete  picture  of 

the  cathode  layer,  the  region  where  electron  avalanching  multiplies  the  posi¬ 
tive  ion  driven  cathode  electron  emission  up  to  the  level  needed  to  carry  the 
discharge  current,  which,  at  the  cathode,  is  primarily  an  ion  current.  What 
the  standard  theory  shows  is  that,  for  a  given  gas  and  cathode  material,  the 
cathode  layer  voltage  drop  is  essentially  fixed,  independent  of  pressure  (more 
precisely,  of  gas  density).  The  sheath  voltage  drop  is  just  that  needed  to 
give  the  required  lO-to-20  fold  electron  multiplication  in  the  sheath.  In 
nitrogen  this  sheath  voltage  is  about  500  V. 

The  electric  field  in  the  sheath  is  directly  proportional  to  density  and 
the  sheath  thickness  is  inversely  proportional  to  density;  i.e.,  there  is  a 
constant  voltage  drop  per  electron  mean  free  path  and  a  constant  number  of 
mean  free  paths  within  the  sheath.  Since  the  current  at  the  cathode  is  mostly 
ion  current,  there  is  a  net  positive  charge  density  within  the  sheath  and, 
since  this  must  satisfy  Poisson's  equation  for  the  voltage  drop,  the  sheath 
has  a  characteristic  current  density,  which  varies  inversely  as  the  square  of 
the  sheath  thickness.  So  the  current  is  proportional  to  the  square  of  the  gas 
density.  At  high  pressures,  this  current  density  is  greater  than  the  dis¬ 
charge  current,  so  the  sheath  breaks  up  into  cathode  spots. 

The  foregoing  is  all  quite  standard  analysis  but  worth  reviewing  because 
it  is  crucial  to  the  problem  of  streamer  initiation.  The  formation  of  cathode 
spots  gives  an  unavoidable  nonuniformity  that  cannot  be  eliminated  by  contour¬ 
ing  the  electrodes.  In  some  of  the  framing  camera  photographs,  the  spots  are 
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clearly  visible  and  it  appears  likely  that  they  are  the  source  of  streamers 
that  are  then  seen  to  develop.  In  a  switch  discharge,  the  current  density 
will  be  higher,  and  one  might  think  that  this  would  give  more  complete  cover¬ 
age  of  the  cathode  by  the  sheath.  But  the  scaling  analysis  in  Section  2.0 
showed  that,  to  make  a  practical  switch,  one  wants  to  work  at  higher  pressures 
(10  atmospheres  or  more).  Since  the  cathode  layer  current  density  scales  as 
gas  density  squared,  the  concentration  of  the  cathode  current  into  spots  can, 
unfortunately,  be  expected  to  be  worse,  not  better,  than  in  laser  discharges. 

The  other  difference  between  low-  and  high-pressure  discharges  is  that 
heating  of  the  cathode  layer  is  a  much  stronger  effect  at  high  pressure. 

Since  the  current  density  increases  as  gas  density  squared,  and  the  electric 
field  is  proportional  to  the  gas  density,  the  dissipation  in  the  sheath  in¬ 
creases  as  the  cube  of  the  room  temperature  pressure  of  the  discharge.  In 
high  pressure  discharges,  the  cathode  layer  rapidly  heats  and  expands.  The 
gas  density,  and  hence  current  density,  goes  down,  so  the  cathode  spots  can 
expand  to  cover  the  whole  cathode.  Indeed,  some  of  the  experiments  showed 
evidence  of  this.  In  some  shadowgraph  framing  camera  sequences,  one  can  see 
an  expanding  cathode  layer  which  appears  to  cover  the  whole  cathode.  But 
streamers  still  formed  and  the  discharge  still  arced. 

The  reason  that  streamers  still  formed  is  indicated  by  another  shadow¬ 
graph  sequence  (shown  in  Ref.  8)  where  a  sinusoidal  perturbation  of  the  cath¬ 
ode  front  is  visible.  Apparently  the  cathode  layer  expansion  had  become 
unstable,  and  an  analysis  of  the  sheath  expansion  explains  why  this  happens. 
Since  our  analysis  of  sheath  expansion  is  given  in  Ref.  8,  it  need  only  be 
summarized  here.  The  essential  point  is  that,  as  the  layer  expands,  its  im¬ 
pedance  drops.  When  the  cathode  spots  have  merged  to  cover  the  whole  cathode, 
the  sheath  current  density  is  just  that  of  the  bulk  discharge.  Then,  as  the 
layer  expands,  the  sheath  electric  field  decreases,  producing  an  increasing 
cathode  layer  conductivity.  When  the  sheath  conductivity  exceeds  that  of  the 
bulk  discharge,  which  typically  happens  at  a  thickness  of  about  1  ram,  the 
expansion  becomes  unstable.  When  the  cathode  layer  is  a  better  conductor  than 
the  bulk  discharge,  a  bump  on  the  cathode  layer  will  draw  an  increased  cur¬ 
rent,  causing  further  heating  and  more  rapid  expansion  at  that  location.  This 
generates  a  streamer.  In  fact,  this  mechanism  is  quite  similar  to  the  process 
responsible  for  streamer  growth. 
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There  Is  no  evident  way  to  stabilize  the  cathode  layer  expansion. 

Even  If  streamer  growth  of  cathode  spots  could  be  prevented,  which  In  a  high 
pressure  switch  discharge  seems  unlikely,  one  must  still  assume  that  streamers 
will  be  Initiated  at  the  cathode.  The  one  thing  that  could  be  done  Is  to  use 
a  porous  or  mesh  cathode  and  flow  gas  Into  the  cathode:  I.e.,  suck  the  cathode 
layer  and  the  associated  small  streamers  Into  the  cathode  surface  fast  enough 
to  prevent  streamer  growrth  Into  the  discharge.  The  needed  flow  Is  just  the 
velocity  at  which  a  streamer  grows  early  In  Its  development.  This  solution, 
however,  would  make  the  switch  apparatus  much  more  complicated  and  expensive. 
The  velocity  of  a  small,  beginning  streamer,  while  less  than  that  of  larger 
ones,  is  still  considerable  -  of  the  order  of  0.1  time  the  sound  speed  In  the 
gas.  This  solution  should  be  kept  In  mind,  but  only  as  a  last  resort.  If 
possible,  the  switch  should  be  designed  to  prevent  arcing  by  keeping  the  pulse 
length  shorter  than  the  time  needed  for  streamers  to  propagate  across  the  dis¬ 
charge. 

3.3  THEORIES  OF  STREAMER  GROWTH 

The  formation  of  arcs  In  e-beam  sustained  discharges  has  not  been  stud¬ 
ied  nearly  as  Intensively  as  arc  formation  In  self-sustaining,  overvoltaged 
discharges,  such  as  preionized  TEA  lasers,  spark  gaps,  lightning  strokes, 
etc.  But,  as  reported  in  Refs.  5  and  6  and  the  papers  cited  therein,  some 
calculations  have  also  been  made  In  an  attempt  to  understand  the  arcing  of 
nonself-sustained  discharges.  However,  comparison  with  the  observations  dis¬ 
cussed  above  does  not  support  these  previous  analyses. 

The  predictions  of  arcing  time  are  In  some  cases  of  the  right  order,  but 
the  observed  dependence  upon  discharge  voltage  and  current  Is  not  predicted. 
The  calculations  often  give  a  characteristic  streamer  velocity,  analogous  to 
the  overvoltaged  case,  or  else  predict  a  collapsing  current  throughout  the 
discharge,  not  the  growth  of  filaments  from  the  electrodes  that  Is  seen  In 
actual  devices. 

A  conmon  procedure  In  stability  analyses  Is  to  linearize  the  fundamental 
equations  by  assuming  small  perturbations  of  electron  density,  gas  density, 
etc.,  within  a  uniform  discharge.  Then,  In  the  same  way  that  one  analyzes 
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sound  waves  in  a  gas,  or  the  various  waves  and  instabilities  in  a  fully  ion¬ 
ized  plasma,  one  can  obtain  a  set  of  linear  equations  for  the  perturbation  and 
study  its  stability.  If  small  disturbances  are  found  to  grow  in  amplitude, 
the  discharge  is  predicted  to  be  unstable. 

The  problem  with  this  technique  is  that  it  assumes  that  the  perturba¬ 
tions  are  initially  of  small  amplitude,  so  that  terms  quadratic  in  the  varia¬ 
tions  in  densities,  field,  etc.,  can  be  neglected.  But  the  experiments  on 
streamer  growth  do  not  support  such  a  picture.  The  photographs  show  a  large 
amplitude  disturbance  that  grows  in  size.  The  variation  in  gas  density,  cur¬ 
rent  density,  etc.,  within  the  streamer  is  not  a  small  perturbation,  even 
early  in  the  process.  So  amplitude  linearization  is  not  an  appropriate  way  to 
describe  this  phenomenon,  even  approximately. 

Existing  analyses  also  tend  to  consider:  increased  local  gas  heating 
caused  by  current  concentration,  and  increased  local  electron  multiplication 
caused  by  E  field  enhancement  as  two  different  instability  modes  to  be  dis¬ 
cussed  separately.  But,  as  explained  above,  both  the  visual  appearance  of 
streamers  and  the  functional  dependence  of  their  growth  rate  upon  the  dis¬ 
charge  parameters  strongly  suggest  that  gas  heating  and  electron  avalanching 
are  both  essential  to  streamer  propagation.  It  appears  that  the  two  processes 
must  be  considered  together  to  construct  a  meaningful  theory. 

Both  effects  were  considered  in  an  earlier  theoretical  model  that  we 
constructed  to  describe  streamer  growth  in  laser  discharges.  Since  the  work 
done  under  the  present  program  is  a  generalization  of  that  earlier  analysis, 
we  will  first  summarize  the  simple  model.  (It  is  discussed  more  fully  in 
Refs.  7  and  8.)  Jhe  basic  idea  is  that  a  streamer  is  a  cylinder  of  hot  and 
consequently  conductive  gas  that  protrudes  into  the  discharge.  Near  its  head 
is  both  a  strong  field  and  a  high  current  density,  just  as  would  occur  near 
the  tip  of  a  wire  protruding  from  the  cathode.  Both  field  and  current  act  to 
raise  the  conductivity  of  the  gas  ahead  of  the  streamer — the  field  by  ioniza¬ 
tion  and  the  current  by  heating. 

To  estimate  the  rate  of  growth  of  a  streamer  one  must  compare  the  energy 
input  needed  to  make  the  gas  ahead  of  the  streamer  highly  conducting  with  the 
rate  of  heating  caused  by  the  concentrated  current  ahead  of  the  streamer.  For 
nitrogen,  in  which  most  of  the  experiments  were  done,  the  Saha  equation  pre¬ 
dicts  ionization  and,  hence,  high  conductivity,  at  ~  6200"K.  At  constant 
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pressure,  the  energy  needed  to  heat  a  constant  volume  of  nitrogen,  from  which 

3 

the  gas  is  permitted  to  expand  to  this  temperature,  is  c  ~  1  J/cm  atm. 

(This  may  be  an  underestimate;  in  nitrogen,  a  signifiant  fraction  of  the 
energy  may  be  taken  up  by  vibrational  excitation,  whfch  is  stored  as  internal 
energy.) 

Around  the  tip  of  the  streamer  there  is  a  halo,  of  radius  a,  where  the 
conductivity  is  high  due  to  nonthermal  ionization  by  the  strong  electric 
field.  In  our  original  model  we  assumed  that  most  of  the  heating  therefore 
occurs  at  the  halo  edge.  Modeling  this  by  a  conductivity  that  varies  as 

(Ic;  -1) 

o  =  e  (12) 


where  r  is  the  distance  from  the  tip  of  the  hot  gas  column,  one  can  integrate 
2 

j  /o  over  the  halo  from  the  radius  of  the  column,  b,  to  the  radius  of  the 
halo,  a,  to  obtain  the  total  dissipated  power  in  the  halo. 
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Since  the  ratio  of  conductivities  in  the  core  at  6200*K  and  in  the  background 
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gas  under  e-beam  ionization  is  ~  10  and  a/b  w  3,  the  value  of  k  which 
matches  the  boundary  condition  on  the  conductivity  is  k  =  3.45.  Thus,  the 
second  term  In  the  above  equation  is  completely  negligible  and  the  mean  power 
density  within  the  halo  is 
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q  =  - 5 - 7^  W/cm 
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From  this  power  density,  one  could  compute  the  time  required  to  heat  the 
halo  to  the  temperature  of  the  core  of  the  streamer.  However,  a  heating  time 
calculated  In  this  way  would  be  too  long  because  only  the  center  need  be 
heated  for  the  core  to  grow.  One  can  show  that,  to  be  consistent  with  the 
assumed  radial  dependence  of  the  conductivity,  the  ratio  of  the  volume  to  be 
heated  to  the  final  temperature  should  be  reduced  by  a  Factor  X  >  0.3. 

These  results  Imply  an  Instantaneous  velocity  for  propagation  of  the 
streamer.  It  Is  the  radius  of  the  sphere  a.  divided  by  the  time  required  to 
heat  the  gas  In  the  halo  Interior  to  the  core  temperture.  t  >  Xc/q.  Thus  we 
have 
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where  I  Is  the  total  current  through  the  streamer  and  c  Is  the  streamer 
length.  This  current  can  be  estimated  by  assuming  that  the  streamer  Is  essen¬ 
tially  at  cathode  potential,  that  the  current  Is  quasi-static  (I.e..  diver¬ 
genceless).  and  that  the  conductivity  outside  the  halo  Is  uniform.  Then  the 
potential  beyond  the  streamer  obeys  Laplace's  equation. 

▼  •j  ^  V  -  0 


Approximating  V  by  the  first  few  Legendre  polynomials,  matching  the  boundary 
conditions  on  V  at  points  on  the  halo  and  the  anode,  assuming  J  »  at  the 
halo  edge,  and  Integrating  over  the  cap  sphere  gives  a  total  current 
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where  E^  Is  the  electric  field  far  from  the  streamer.  Using  this  In  the 
expression  for  the  velocity  gives  a  propagation  velocity  for  the  streamer 
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The  model  predicts  that  c/a  should  remain  roughly  constant  as  the 
streamer  grows  and  this  Is  confirmed  by  the  framing  camera  photographs.  Put¬ 
ting  the  result  In  a  form  more  appropriate  for  such  a  scaling,  and  using  the 
above  estimated  values  of  k,  and  c.  we  have 

dc 


where 


Hence  the  model  predicts  that  streamers  should  grow  exponentially.  This  pre¬ 
diction  Is  In  good  agreement  with  the  observations.  The  predicted  growth 

2 

rate,  y.  Is  primarily  dependent  upon  the  specific  power  loading,  P^. 

This  Is  also  consistent  with  the  data  at  the  relatively  high  electric  fields 
at  which  laser  discharges  run. 

However,  as  discussed  above,  data  taken  at  lower  fields  more  typical  of 
switch  discharges  show  a  different  parametric  variation.  The  most  probable 
reason  that  the  foregoing  model  falls  there  Is  that  heating  ahead  of  the  halo 
was  neglected,  and  such  heating  Is  more  significant  at  lower  fields,  where  the 

halo  will  be  relatively  small.  To  estimate  the  effect  of  heating  gas  outside 

the  halo,  one  can  consider  the  other  extreme  and  assume  that  streamer  growth 
Is  caused  by  this  alone.  Just  beyond  the  halo,  the  electric  field  Is  E^, 

the  critical  breakdown  field,  and  the  conductivity  Is  o  .  So  the  power  1n- 

2  ° 
put  Is  OqEj.#  which  will  heat  the  gas  to  the  conductive  temperature  In  a 

time  t  *  e/o^E^.  Since  the  scale  length  of  this  field  Is  a,  the  halo 

radius,  such  heating  alone  would  give  a  velocity  v  -  a/t.  Here  a,  the  halo 

radius,  can  be  estimated  from  the  Legendre  polynomial  expression  for  the 

E-fleld  mentioned  earlier.  The  result  Is: 


c 
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This  gives  a  speed 
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(20) 


dt 
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So  one  again  predicts  exponential  growth,  but  now  the  growth  rate  Is* primarily 
dependent  upon  o^E^  =  the  current  density  In  the  uniform  background 
discharge.  This  Is  an  extreme  conclusion,  since  It  Ignores  all  heating  within 
the  halo,  but  It  Is  worth  noting  that  the  data  In  Figure  3  shown  earlier  are 
not  really  Inconsistent  with  a  growth  rate  that  Is  primarily  current-dependent 
at  low  discharge  voltages. 

Clearly,  what  Is  needed  Is  a  more  complete  description  that  treats  the 
whole  region  around  the  end  of  a  streamer  consistently.  Before  proceeding  to 
that,  however,  we  would  like  to  review  one  other  simple  estimate  that  Is  use¬ 
ful  as  a  guide  to  the  analysis.  As  noted  earlier,  the  electric  field  within 
the  halo  must  be  nearly  equal  to  the  breakdown  threshold  level.  Since  this 

field  limit  Is  proportional  to  gas  density,  n,  and  the  current  density  must 
2 

vary  as  1/r  If  one  assumes  approximate  spherical  symmetry,  the  power  depos¬ 
ition  In  the  halo  must  have  the  radial  dependence. 
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(22) 


where  r  Is  the  distance  from  the  streamer  tip  and  n  and  E  are  the  gas 

0  CO 

density  and  breakdown  field  of  the  background  discharge. 


Matching  this  to  the  boundary  condition  at  the  halo  edge  gives 


^  =  **0  'cS 
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For  nitrogen,  an  analysis  of  gas  heating  shows  that. 


e  *  c. 


0.63 


Differentiating  this  gives. 
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dT. 


We  also  know  that  n  »  n  T  /T,  so 

0  o 


dT  = 


n^T  dn 
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Combining  these  results  gives 
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Making  the  substitution  dn/dt  «  v  dn/dr,  we  obtain 


.-2.63  dn  K 
dr  “  ^2 


where 


T  0.63  .22 

^1  Vco  * 


OM  n.  T.  V 
o  0 


This  has  the  solution,  for  n  =  n  at  r  =  a. 
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Well  Inside  the  halo,  this  reduces  to 


n  =  K' 


0.61 

r 


(30) 


So  one  expects  the  gas  density  to  vary  roughly  as  the  0.6  power  of  the  dis¬ 
tance  from  the  streamer  tip  In  the  Interior  of  the  halo.  The  point  Is  that, 
from  such  an  analysis,  one  can  deduce  dependences  that  had  merely  to  be  postu¬ 
lated  In  the  earlier  model. 

The  objective  now  Is  to  start  from  the  basic  equations  for  gas  motion, 
electric  field,  etc.,  and  try  to  put  the  whole  picture  on  a  solid  mathematical 
footing.  Since  one  cannot  linearize  by  assuming  small  amplitude  perturba¬ 
tions,  that  might  seem  to  be  a  hopeless  undertaking.  But  the  data  suggest 
that  streamer  growth  has  a  different  kind  of  linearity.  Streamers  were  seen 
to  grow  exponentially,  which  Is  a  characteristic  of  a  linear  system.  We  ex¬ 
pect  that  the  equations  will  be  at  least  linear  In  time. 

The  picture  that  emerges  from  the  simple  calculations  and  the  data  Is 

that  a  streamer  grows  In  size  while  retaining  the  same  form.  Such  a  structure 

can  be  described  by  a  similarity  solution,  a  technique  which  Is  used  In  fluid 
( 28) 

mechanics.  One  assumes  that  the  coordinates  of  the  structure  have  an 
overall  time  dependence.  In  this  case  that  they  vary  as  exp(Yt).  The  Idea  Is 
to  try  to  account  for  all  the  time  variation,  or  at  least  all  the  rapid  varia¬ 
tion,  In  this  simple  way.  The  framing  camera  photographs  strongly  suggest 
that  such  a  solution  to  the  equations  should  exist. 


4.0  MATHEMATICAL  MODEL  OF  A  STREAMER 


4.1  ASSUMPTIONS  OF  THE  MODEL 

As  In  our  earlier  work,  a  streamer  In  the  new  model  consists  of  three 
regions:  the  core,  the  halo,  and  the  background  discharge.  These  three  re¬ 
gions  are  distinguished  by  the  source  of  their  electrical  conductivity.  The 
core,  essentially  an  arc.  Is  a  column  of  hot,  low-density  gas  that  Is  therm¬ 
ally  Ionized  and  highly  conducting.  In  this  region  there  Is  a  high  current 
density  and  a  weak  electric  field.  Around  the  tip  of  the  core  Is  the  streamer 
halo,  where  the  electric  field  Is  strong  enough  to  produce  nonthermal,  elec¬ 
tron  avalanche  Ionization  of  the  gas.  Outside  the  halo  Is  the  background  re¬ 
gion  where  the  conductivity  Is  caused  by  e-beam  Ionization.  The  background, 
which  has  a  lower  conductivity  than  the  core  or  halo,  acts  as  a  resistive 
ballast,  limiting  the  current  that  can  flow  through  the  streamer. 

As  sketched  In  Figure  6,  the  core  of  the  streamer  Is  modeled  as  a  cylin¬ 
der  of  length  c  and  radius  b.  The  radius  may  vary  along  the  length,  since 
heat  conduction  will  cause  the  core  to  expand  with  time,  but  this  variation  Is 
so  gradual  that  the  dynamics  are  still  those  of  a  cylindrical  structure. 
Moreover,  all  that  Is  actually  needed  In  the  calculation  Is  the  radius  near 
the  tip  of  the  column. 

The  halo  Is  modeled  as  a  region  of  radius  a  around  the  tip  of  the  core. 
As  shown  by  our  earlier  photographic  studies,  this  region  Is  not  exactly 
spherical,  but  Is  nearly  so,  and  all  that  Is  needed  for  the  calculations  Is 
the  structure  of  the  region  ahead  of  the  streamer,  where  gas  heating  causes 
streamer  growth. 

Some  distinctions  In  terminology  should  be  noted.  The  halo  discussed  In 
the  theory  Is  the  region  of  nonthermal  Ionization  by  the  electric  field.  The 
glowing  ball  seen  In  the  photographs  Is,  presumably,  a  region  of  nonthermal 
luminosity.  Since  an  electric  field  too  weak  to  cause  electron  avalanching 
could  still  accelerate  electrons  to  energies  sufficient  to  excite  optical 
transitions,  the  visible  ball  may  be  somewhat  larger  than  the  Ionization 
halo.  But  the  fact  that  the  ball  appears  spherical  Is  certainly  a  strong  In¬ 
dication  that  this  region  Is  roughly  spherical. 


Note  also  that  the  background  discharge,  as  defined  here,  is  not  neces¬ 
sarily  the  same  as  the  unperturbed  discharge  far  from  the  streamer.  The  back¬ 
ground  Is  simply  the  region  outside  the  halo,  where  Ionization  Is  caused  by 
the  e-beam.  The  gas  In  this  region  ahead  of  the  streamer  may  be  heated  enough 
by  the  streamer  current  to  significantly  change  Its  density  and  conductivity. 

In  the  following  analysis,  the  subscripts  c,  H.  and  B  are  used  to  denote 
core,  halo,  and  background  regions.  Thus  T^  and  are  the  temperature 
and  electrical  conductivity  of  the  core,  o  (|E|  ,  T,  P)  Is  the  conductivity 
of  the  halo,  etc.  The  properties  of  the  unperturbed  discharge  far  from  the 
streamer  are  denoted  by  the  subscript  naught.  Thus  p^,  E^,  and  are 
the  gas  density,  electric  field,  and  current  density  of  the  unperturbed  dis¬ 
charge.  These  quantities  may  be  time  varying,  but  they  are  not  dependent  on 
position.  It  Is  assumed  that  the  whole  pattern  Is  cylindrical ly  symmetric 
about  the  axis  of  the  streamer.  For  pulsed,  uniform  discharges,  this  Is  well 
borne  out  by  the  data. 

It  Is  also  assumed  that  the  local  properties  are  Isotropic,  that  the 
pressure,  p,  electrical  conductivity,  o,  and  thermal  conductivity,  K,  are 
scalars,  not  tensors.  Consistent  with  the  use  of  a  scalar  pressure,  the  vis¬ 
cosity  and  the  resultant  viscous  dissipation  of  energy  are  neglected.  Be¬ 
cause,  as  noted  earlier,  the  data  do  not  show  any  evidence  of  shock  waves  or 
other  sonic  effects  as  the  streamer  velocity  exceeds  the  sound  speed  of  the 
background  gas,  the  gas  flow  away  from  the  streamer  Is  assumed  to  be  sub¬ 
sonic.  Consistent  with  subsonic  flow,  the  pressure  variations  are  expected  to 
be  relatively  small,  and  the  kinetic  energy  of  flow  of  the  gas  Is  expected  to 
be  much  less  than  the  changes  In  Internal  energy  caused  by  heating  by  the 
current. 

Finally,  the  gas  kinetic  processes  -  Ionization,  free  electron  heating, 
etc.,  are  assumed  to  be  rapid  In  comparison  with  the  streamer  growth  rate. 
Hence  the  electrical  conductivity,  for  example.  Is  considered  to  be  a  function 
of  the  Instantaneous  gas  temperature,  density,  and  electric  field,  but  not  to 
depend  upon  the  history  of  the  gas.  This  omits  such  effects  as  E-fleld  over¬ 
shoot  when  the  Ionization  lags  the  current  Increase,  or  non-thermal  storage  of 
energy  In  vibrational  modes,  and  subsequent  abrupt  heating  when  the  gas  tem¬ 
perature  rises  to  the  point  where  vibrations  relax  rapidly.  Such  effects 
could  be  included  through  additional  variables,  such  as  a  vibration  tempera¬ 
ture,  which  could  be  different  from  t^e  gas  temperature,  or  by  allowing  rate 


dependence,  e.g.,  a  conductivity  dependent  upon  dE/dt  as  well  as  E,  but  this 
was  considered  to  be  an  impractical  complication.  The  rapid  turn-off  needed 
In  a  switch  Implies  that  kinetic  rates  will  be  rapid  -  submicrosecond  -  while 
the  needed  conduction  time  implies  that  streamer  growth  rates  must  be  slower  - 
tens  of  microseconds  -  in  systems  of  practical  interest. 

4.2  BASIC  EQUATIONS 
4.2.1  Gas  Dynamics 


Since  the  system  is  cylindrical 1y  symmetric,  it  can  be  described  in 
terms  of  three  independent  variables,  r,  z,  and  t,  with  r  and  z  defined  as 
sketched  in  Figure  6.  Streamers  are  assumed  to  originate  at  an  electrode  and 
propagate  into  a  uniform  discharge.  Because  the  opposing  electrode  affects 
the  streamer  only  when  it  has  nearly  bridged  the  gap,  when  only  a  brief  time 
remains  before  arcing,  most  of  the  streamer  growth  time  can  be  explained  with¬ 
out  reference  to  an  opposing  electrode.  Hence  the  discharge  will  be  consid¬ 
ered  to  be  semi-infinite,  with  the  boundary  condition  that  the  current, 
electric  field,  etc.,  must  approach  the  uniform,  unperturbed  values  at  large 
distances  from  the  streamer. 

The  success  of  our  earlier  theoretical  estimates  in  predicting  the  right 
magnitude,  and,  at  high  electric  fields,  the  right  scaling  for  the  growth  rate 
of  a  streamer,  is  persuasive  evidence  for  the  assumption  that  streamer  growth 
is  due  to  gas  heating.  Because  the  data  also  show  no  evidence  of  shock  waves, 
one  expects  that  the  pressure  around  a  streamer  will  be  nearly  uniform.  For 
these  reasons,  it  is  convenient  to  describe  the  state  of  the  gas  in  terms  of 
temperature,  T,  and  pressure,  p,  and  to  describe  the  internal  energy— heating, 
dissociation,  etc.— by  the  internal  enthalpy  per  unit  mass,  h  =  e  +  p/p,  where 
e  is  the  internal  energy.  The  quantity,  h,  is  assumed  to  be  known  as  a  func¬ 
tion  of  T  and  p  for  the  discharge  gas  mixture.  The  other  relevant  gas  proper¬ 
ties,  the  density,  p.  and  thermal  conductivity,  K,  are  also  assumed  to  be 
known  or  calculable  in  terms  of  T  and  p. 

)r  addition  to  T  and  p  (and  p,  h,  and  K),  there  are  six  other  variables, 
the  two  components  of  velocity,  v^  and  v^,  current,  and  j^,  and 
electric  field,  E^  and  E^.  Eight  equations  are  therefore  needed  to 
describe  the  dynamics  of  the  system. 

One  knows  that  the  mass  and  momentum  of  the  gas  are  conserved. 


+  v-(pv)  *  0 
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Here  0/Dt  denotes  the  comovlng  derivative,  a/at  +  y. r ,  and  (32)  includes  the 
fact  that  there  are  no  forces  on  the  gas,  and  that  the  viscous  drag  is  neglig¬ 
ible. 

The  effect  of  gas  heating  by  the  current  is  described  by  the  enthalpy 
equation. 


^  Ot  "  ft  *1-1  (33) 

Here  viscous  dissipation  has  been  neglected  and  the  equation  has  been  written 
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in  terns  of  the  total  enthalpy,  H  =  v  /2  +  h.  Since  the  gas  flow  is  sub¬ 
sonic  and  the  gas  heating  Increases  h  manyfold,  one  expects  that  the  differ¬ 
ence  between  H  and  h  will  probably  be  unimportant. 

4.2.2  Electrodynamics 

Because  there  are  no  significant  magnetic  fields,  one  has,  from 
Faraday's  law. 


▼  X  E  »  0 


(34) 


Poisson's  equation  and  change  conservation  together  imply  that 

V,  [1*  -  0 

Because  the  charge  transported  by  the  current  in  these  discharges  is 
many  orders  of  magnitude  greater  than  the  charge  needed  to  produce  the  elec¬ 
tric  field,  the  second  term  is  quite  negligible,  and  one  has,  to  good  approx¬ 
imation. 
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since  the  momentum  equation,  (£q.  32)  has  two  components.  Eqs.  (31-35) 
provide  six  independent  equations.  The  remaining  two  that  are  needed  are,  of 
course,  the  Ohm's  law  relation  between  i  and  E.  As  discussed  already,  this 
dependence  is  different  for  the  different  regions  of  the  streamer.  In  all 
three  regions,  however,  the  conductivity  is  scalar,  so  one  always  has 


(35) 


All  that  is  still  needed  is  the  dependence  of  |  j  |  upon  |  e|. 

4.2.3  Conductivity 
a.  Background 

In  the  background  region,  the  conductivity,  o.  -  n  n  e  is  propor- 

D  e  e 

tional  to  the  electron  density,  n^,  which  is  determined  by  the  balance 
between  e-beam  ionization  of  the  gas  and  electron  loss  through  recombination 
or  attachment.  The  electron  loss  rates  may  depend  upon  the  electric  field 
(through  the  free  electron  temperature).  The  mobility,  is  also  in  gen¬ 
eral  a  function  of  E/#>.  Thus  the  only  generally  correct  form  is  o-  =  o 

B  B 

(T.p,  E  ),  or,  solving  i  «  o  E  for  E(i), 

Ob  •  <»B  <37a)) 

In  many  cases  of  Interest,  however,  the  mobility  and  the  recombination  or 
attachment  rates  may  be  approximated  as  independent  of  E,  because  the  electron 
temperature  is  clamped  by  the  energy  sinks  provided  by  atomic  or  molecular 
transitions.  If  the  gas  is  well  below  the  ionization  temperature,  these  rates 
may  depend  on  T  and  p  only  through  p.  Then,  since  the  rate  of  production  of 
electrons  by  the  e-beam  is  proportional  to  p,  one  has  the  simpler  dependence, 

gg  =  Og(p)  (37b) 


Of  course,  a  dependence  on  e-beam  current  is  also  implied  here. 
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If  the  electron  loss  Is  dominated  by  attachment,  whose  rate  Is  propor¬ 
tional  to  p,  then.  Instead  of  Eq.  (37c),  one  has  n^  -  n^^,  and  hence 


b.  Halo 

In  the  halo,  the  Ionization  Is  enhanced  by  electron  avalanching.  Since, 
near  the  breakdown  threshold,  this  rate  Is  a  rapidly  Increasing  function  of 
electric  field.  It  Is  Impractical  to  try  to  estimate  n^  from  the  electron 
production  and  loss  rates.  Fortunately,  however.  It  Is  not  necessary  to  do 
that.  Because  the  rate  of  electron  production  Increases  rapidly  as  E  exceeds 
the  breakdown  strength  of  the  gas,  one  expects  that  such  an  Increase  would 
rapidly  (In  a  time  short  compared  to  streamer  growth  times)  Increase  the  con¬ 
ductivity,  which  would.  In  turn,  cause  the  electric  field  to  drop  back  towards 
the  breakdown  field.  (Remember  that  the  streamer  current  Is  ballasted  by  the 
resistance  of  the  background  discharge;  so,  near  the  streamer,  the  system  Is 
effectively  driven  by  a  current  source.)  Hence  we  hypothesize  that.  In  the 
halo,  the  electric  field  Is  clamped  at  the  breakdown  threshold: 


^bo  ^0 
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where  E.  Is  the  breakdown  strength  of  the  unperturbed  discharge.  If  de- 
bo 

sired.  Eg.  (38a)  can  be  expressed  In  terms  of  conductivity,  o^(p,j),  by 
dividing  J  by  E. 


This  expression  for  the  breakdown  field  Ignores  any  gas  temperature 
dependence  that  could  arise  If  the  gas  were  hot  enough  for  a  significant  frac¬ 
tion  of  the  atoms  to  be  In  excited  states,  but  that  Is  expected  to  be  a  small 
correction,  especially  since  the  system  Is  at  nearly  constant  pressure:  as  the 
gas  Is  heated.  Its  density  decreases,  so  subsequent  heating  acts  upon  a 
smaller  mass  of  gas.  Thus,  as  the  gas  approaches  the  Ionization  temperature, 
the  heating  occurs  more  and  more  rapidly.  Once  the  gas  begins  to  Ionize 
thermally,  the  temperature  will  plateau,  as  additional  energy  Is  Invested, 
mostly  In  Ionization. 

There  may,  however,  be  a  second,  lower  temperature  plateau,  probably 
within  the  halo,  where  the  gas  Is  dissociating.  To  Include  such  effects  one 
could  add  an  explicit  temperature  dependence  to  0^.  That  could  be  done 

without  changing  the  formulation  developed  here.  In  general,  the  division  of 

« 

the  streamer  Into  core,  background,  and  halo  Is  Itself  a  simplifying  approx¬ 
imation.  Near  the  background-halo  boundary,  there  must  be  a  region  where 
e-beam  Ionization  and  electron  multiplication  are  both  significant.  Simi¬ 
larly,  near  the  halo-core  boundary  there  must  be  a  region  where  both  avalanch¬ 
ing  and  thermal  Ionization  are  significant.  The  model  could  be  generalized  by 
constructing  an  overall  conductivity  o(T,  p,  j),  which  Includes  all  these  pro¬ 
cesses,  but  that  has  been  left  for  later  work. 

c .  Core 

In  the  core  of  a  streamer,  the  gas  Is  thermally  Ionized,  as  In  an  arc. 
The  Saha  equation  shows  that  this  happens  fairly  abruptly,  as  the  temperature 
reaches  a  threshold.  Hence  we  expect  that  the  core  Is  at  a  well  defined  tem¬ 
perature 
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where  1s  the  Ionization  temperature  of  the  gas.  That  Is  from  an  estimate 
of  the  degree  of  Ionization  one  can  calculate  the  gas  temperature  from  the 
saha  equation.  In  general,  this  depends  upon  pressure,  but,  since  pressure 
variations  are  expected  to  be  small.  It  can,  to  reasonable  accuracy,  be 
evaluated  at  the  ambient  pressure  of  the  discharge. 

The  degree  of  Ionization  within  the  core,  which  can  vary  with  radius.  Is 
difficult  to  estimate.  Hence  the  conductivity  and  electric  field  In  the  core 
are  not  as  easily  calculated  as  In  the  other  regions.  But  since  the  core  Is 
essentially  an  arc,  one  can  make  some  use  of  existing  understanding  of  arc 
discharges  and,  moreover,  the  exact  characteristics  of  this  region  do  not 
strongly  influence  the  conclusions  about  streamer  growth.  The  core  conductiv¬ 
ity  Is  certainly  high  enough  so  that  the  tip  Is  close  to  cathode  potential. 

One  can  simply  take  a  rough  estimate  obtained,  for  example,  from  arc  physics 
for  this  high  conductivity 


O,  =  o,(Tj)  (39) 

It  would  also  be  useful  to  estimate  the  energy  Invested  In  Ionization  and  the 
radius  of  the  core,  since  this  determines  the  current  density.  But,  because 
the  gas  even  In  the  core  Is  only  weakly  Ionized  and  because  the  final  heating 
near  the  core  tip  is  done  on  gas  of  much  lower  density  from  the  background 
discharge  gas,  these  parameters  also  have  only  a  small  effect  on  streamer 
dynamics. 

This  completes  the  formulation  of  the  problem.  To  solve  these  equa¬ 
tions,  It  is  not  actually  necessary  to  assume  a  structure  with  a  cylindrical 
core  and  a  ball-shaped  halo  in  a  background  discharge.  One  could  regard 
Eqs.  (37-39)  as  defining  a  conductivity  o(T,p,j)  according  to  the  prescrip¬ 
tion:  If  T  >  Tj,  Eq.  (39)  applies.  For  colder  gas,  Eq.  (37)  should  be 
used,  unless  the  predicted  electric  field  would  exceed  the  breakdown  strength 
of  the  gas  at  the  local  density.  In  which  case,  Eq.  (38)  applies.  If  this 
model  Is  correct,  a  structure  like  that  seen  In  the  framing  camera  photographs 
should  emerge  from  a  general  solution  of  these  equations. 

Even  more  Ideally,  one  might  prefer  to  solve  the  equations  numerically 
In  three  dimensions,  to  confirm  the  expectation  of  cylindrical  symmetry.  This 
could.  In  principle,  be  done.  The  addition  of  a  third  spatial  dimension  adds 
three  variables,  the  third  components  of  y,  E,  and  1,  but  It  also  adds  three 
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equations,  the  third  conponent  of  Ohm's  law  and  two  more  components  of 
Faraday's  law.  We  have  thus  succeeded  In  formulating  the  basic  physics  of  the 
streamer  Instability  as  a  well  posed  mathematical  problem,  which  had  not  been 
done  prior  to  this  program. 

4.3  DESCRIPTION  IN  CYLINDRICAL  COORDINATES 

Restating  the  equations  In  cylindrical  coordinates,  with  no  angular 
dependence,  and  only  r  and  z  components  of  the  vectors,  one  has  the  conserva¬ 
tion  of  mass. 


«  *■  r  f;  «  <"z)  ■  “ 

(40) 

where  p  =  p  (T,p),  the  conservation  of  momentum. 
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(41) 

av  av  av  , 

v^  ^  ♦  v_  ^  *  0 

dv  dP  z  dz  p  9z 

(42) 

The  enthalpy  equation 
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where 
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2  *  2  h(T,p) 

Faraday's  law. 


charge  conservation. 


1  a 

?  i7  "  IT  -  0  (45) 

and  Ohm's  law 

=  o(T.p.j)  (46) 

*  «»(T.p.j)  E^  (47) 

where 

j  - 

4.4  CONVERSION  TO  COMOVING.  SPHERICAL  COORDINATES 

As  explained  earlier,  we  propose  to  look  for  a  similarity  solution  to 
these  equations,  a  solution  describing  a  streamer  that  grows  In  size  while 
retaining  the  same  form.  A  convenient  way  to  do  this  Is  first  to  transform  to 
spherical  coordinates  (n,  6)  centered  on  the  tip  of  the  steamer  column,  as 
shown  In  Figure  6.  In  these  terms,  the  cylindrical  coordinates  (r,  z)  are 
given  by 

r  «  n  ^  cos  e  (48a) 

z  -  c  +  n  f  sin  e  (48b) 

Here  c  Is  the  length  of  the  streamer  column,  f  Is  a  scale  factor,  and  c(t)  and 
f(t)  will  be  chosen  to  simplify  the  equations.  Of  course,  one  also  must 
transform  the  vector  components  according  to 

V  *  V  sin  8  +  v^  cos  8  (49a) 

r  n  8 


V  *  V  cos  8  -  V  sin  8  (49b) 

Z  ft  9 

and  similarly  for  the  components  of  J  and  E. 

After  considerable  algebra  (omitted  here),  one  obtains  the  following 
forms.  The  equation  of  mass  conservation  (40)  becomes 
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From  Euler's  Eqs.  (41,42)  one  has  the  two  components 
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The  enthalpy  Eq.  (43)  becomes 
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Faraday's  law  In  the  new  coordinates  Is 

an  n  n  a  O 


(54) 


Charge  conservation  has  the  form 

;  ^  nh  il-  <Je  *>  -  “  <“> 

Ohm's  law  now  has  the  n  and  e  components 

j  =  <»  (56) 

n  n 

and 

Since  the  gas  density  and  thermal  conductivity  are  determined  by  the  tempera¬ 
ture  and  pressure,  p  >  p  (T,p).  K  =  K  (T.p)  and  the  electrical  conductivity, 
as  explained.  Is  also  a  function  of  the  other  variables,  o  -  o(T,p,J),  there 
are  still  eight  dependent  variables,  T,  p  and  the  n  and  e  components  of  v,  J, 
and  E,  which  obey  these  eight  equations. 

From  the  data  and  the  simple  physical  agreements  discussed  above,  we  ex¬ 
pect  that  the  problem  should  have  a  solution  with  an  exponential  time  depend¬ 
ence.  That  Is,  the  equation  should  be  solved  by  defining 


and  assuming,  as  a  trial  form  that 
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(59a) 
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Here,  v  and  v  are  functions  of  n  and  e,  but  not  of  T.  We  postulate 

f|0  90 

that  this  Is  also  true  of  the  remaining  dependent  variables,  I.e.,  that  In 
these  coordinates 


at 


=  0 


(60) 


and,  likewise  for  p,  T,  j  and  E. 

What  one  finds  Is  that  this  form  Is  a  solution  If  the  thermal  conduction 
can  be  neglected. 


K  =  0  (61) 

This  Is  a  reasonable  approximation  beyond  the  streamer  core.  The  arc-like 
core  has  a  diameter  that  Is  determined  by  heat  conduction,  so  that  will  have 
to  be  analyzed  separately  and  probably,  approximately.  Analysis  of  heat  flow 
In  an  arc  Is  complicated,  even  for  the  standard,  stationary  arc,  and  the  pres¬ 
ent,  time-dependent  problem  may  be  even  more  complicated.  We  note,  however, 
that  arcs  are  known  to  run  at  a  roughly  constant  current  density  over  a  con¬ 
siderable  range  of  total  currents.  That  Is,  as  the  current  Is  Increased,  the 
arc  cross  section  Increases  proportionately.  (See,  for  example.  Section  9.7 
of  the  book  by  Cobine,  Ref.  26).  In  a  streamer,  whose  current  Is  propor¬ 
tional  to  the  square  of  Its  length,  such  a  variation  would  make  the  core  diam¬ 
eter  proportional  to  streamer  length;  i.e..  It  would  cause  this  part  of  the 
structure  to  obey  a  similarity  solution,  too. 

In  any  case,  outside  the  core,  the  thermal  gradients  are  very  much 
reduced,  because  the  scale  Is  larger  and  the  temperature  differences  are 
smaller.  So  we  can  assume  that.  In  the  halo  and  the  background  the  gas  heat¬ 
ing  Is  balanced  by  absorption;  I.e.,  an  Increasing  gas  temperature,  and  by 
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convection,  as  heat  Is  carried  off  by  the  expanding  gas,  but  that  conduction 
Is  a  minor  effect,  which  can  be  neglected. 

Under  this  assumption,  one  obtains  a  simpler  set  of  equations.  Conser¬ 
vation  of  mass.  Eg.  (50),  becomes 

.  Y  sin  o  . 

(Tn  -  T  cos  e)  t  - - - 5| 
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I)  o 

Here  the  time-dependence  has  completely  disappeared,  as  desired. 

The  two  components  of  Euler's  Eq.  (51,  52)  become 

a  V 

’  *00  *  r;r  ‘’.o  -  T  n  -  T  cos  e) 


(v^  ♦  T  c  sin  e)  av  ^ 

^ _ —0^ _  /_JlO  _ 


( 


ae 


'eo* 


tHiSL 


(63) 


and 


av. 


T  \ 


eo 


eo 


an 


(V 


no 


-  Tn  - 


T  c^  cos  0) 


♦  —22 - 0 -  ^  ^  j 

n  '30  no' 


el!!!  36 

p  9B 


(64) 


These  equations  retain  a  time-dependence  In  the  pressure  terms.  At  first 
glance,  this  might  appear  to  Invalidate  the  approach,  since  we  are  assuming 
that  the  density,  temperature,  etc.,  depend  only  upon  n  and  0.  Euler's  equa¬ 
tion  shows  that  this  cannot  be  true  of  pressure  grad1e..ts.  The  reason  Is  that 
the  velocities,  as  functions  of  (n,  0),  have  been  assumed  to  Increase  as  exp 
(yt).  Hence  the  acceleration  of  the  gas  must  have  a  similar  dependence.  So 
the  pressure  gradients  must  also  vary  this  way  to  produce  the  needed  accelera¬ 
tion  of  the  gas  away  from  the  streamer.  Moreover,  ap/3n  is  not  a  gradient. 
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since  the  coordinate  n  Is  Itself  expanding.  The  result  Is  that  this  deriva¬ 
tive  must  Increase  as  exp  (2Yt)  to  conserve  gas  momentum. 

Recall,  however,  that  we  are  assuming  that  the  system  Is  nearly 
Isobaric,  because  the  flows  are  subsonic. 

P  =  Pq  +  P' 

where 

p*  «  p 

0 

and  p^  Is  just  the  constant,  uniform  background  pressure  of  the  discharge. 

This  Is  very  different  from  the  behavior  of  the  density,  temperature,  and 
Internal  enthalpy,  which  vary  by  large  factors  within  a  streamer.  Therefore, 
In  figuring  such  quantities  as  p  (T.p)  one  can  to  this  approximation  use  p^ 
and  consider  p  to  be  a  function  only  of  T.  Equations  (63)  and  (64)  then  spec¬ 
ify  what  the  small  variation  In  pressure  must  be.  If  It  turns  out  to  be  com¬ 
parable  to  p  ,  which,  of  course.  It  will  eventually,  since  It  grows  as 
0 

exp  (2Yt).  this  calculation  ceases  to  be  valid.  That  1s  just  the  statement 
that  eventually  the  streamer,  growing  exponentially,  must  become  supersonic. 

At  that  time,  the  streamer  growth  will  apparently  cease  to  be  a  simple  expon¬ 
ential.  But  the  observations  Imply  that  this  Is  unlikely  to  occur  In  systems 
of  practical  Interest.  Moreover,  even  If  It  did  occur  on  the  last  e-fo1d1ng, 
most  of  the  streamer  growth  time  Is  spent  In  the  earlier  stage,  where  the 
streamer  Is  small  and  moving  slowly.  So  the  present  calculation  should  still 
give  a  fairly  accurate  prediction  of  the  discharge  arcing  time. 

Regarding  Eqs.  (63)  and  (64)  as  specifying  pressure  variations  that  are 
so  small  they  are  unimportant  eliminates  one  variable,  the  changing  pressure, 
and  one  equation.  The  two  equations  still  leave  one,  which  can  be  obtained  by 
taking  the  cross  derivatives  and  differencing  to  eliminate  p.  This  resulting 
equation,  while  quite  complicated,  and  of  second  order.  Is  time-independent  as 
desired. 

The  enthalpy  Eq.  (S3)  becomes 
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Here  heat  flow  and  the  enthalpy  changes  due  to  pressure  variations  have  been 
omitted  because  they  are  small  effects,  and  the  result  Is  an  equation  with  no 
explicit  time-dependence.  Note  also  that  If  E/p  Is  constant,  which  we  believe 
It  to  be  within  the  halo,  where  the  E  field  Is  at  the  avalanche  threshold,  the 
right-hand  side  of  this  equation  Is  just  proportional  to  j,  the  current 
density. 

The  remaining  Eqs.  (54-57),  namely  Faraday's  law,  charge  conservation, 
and  Ohm's  law,  are  unchanged  because  they  were  already  time  Independent. 

We  have  now  obtained  a  set  of  seven  time  Independent  equations  for  the 
seven  remaining  variables.  This  confirms  our  original  supposition  that  a 
similarity  solution  should  exist.  The  next  step  Is  to  simplify  this  set  to 
derive  a  prediction  for  y,  the  growth  rate  of  the  streamer. 

4.5  THE  DYNAMICS  DIRECTLY  AHEAD  OF  THE  STREAMER 

The  remaining  equations,  as  they  stand,  are  still  too  complicated  to 
solve,  at  least  within  the  scope  of  the  present  program.  Thus,  although  the 
most  desirable  result  would  be  a  solution  of  the  whole  set  predicting  the 
shape  of  the  streamer,  the  current  distribution  within  and  around  the 
streamer,  etc.,  we  will  not  be  able  to  carry  the  calculation  to  that  point. 
However,  what  Is  really  needed  for  applications  Is  a  calculation  of  the  growth 
rate  y  and,  for  that  purpose,  It  Is  sufficient  to  consider  only  the  region 
along  the  axis  ahead  of  the  streamer. 

Restricting  the  problem  to  the  line  e  =  0  ahead  of  the  streamer  Is  a 
considerable  simplification.  There,  the  field,  current,  and  gas  flow  must,  by 
symmetry,  be  radial.  E  ,  J  ,  V  and  their  n  derivatives  all  vanish, 
which  gives  a  much  simpler  set  of  equations.  Conservation  of  mass,  Eq.  (62), 
becomes 
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The  n  component  of  Euler's  Eq.  (63)  becomes 
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The  e  component  of  Euler's  Eq.  (64)  is  satisfied  identically,  because 
a11  the  terms  vanish  along  6  =  0.  The  enthalpy  Eq.  (65)  becomes 

dH  1 

?n  <\o  -  ^  Jn  ®  ° 

Faraday's  law  becomes 

aE 

^  =  0  at  0  =  0  (69) 


Charge  conservation  becomes 

i  ^  2 “  •» » ' »  <"» 

and  the  remaining  single  component  of  Ohm's  law  is 

j  =oE  (71) 

n  n 

As  they  stand,  these  six  equations  (actually  five,  since  the  remaining 
Euler  equation  just  gives  the  pressure  variation,  which  is  a  higher-order 
quantity)  are  not  a  complete  set,  as  a  check  of  the  number  of  variables  easily 
shows.  Indeed,  one  would  not  expect  to  be  able  to  solve  a  set  of  partial  dif¬ 
ferential  equations  by  restricting  to  a  line  in  this  fashion.  However,  this 
set  is  sufficient  for  an  estimate  of  the  growth  rate  y  if  one  knows,  or  is 
willing  to  assume,  enough  about  the  form  of  the  solution. 


The  needed  additional  assumption  is  that,  ahead  of  the  streamer,  the 
system  is  nearly  spherically  symmetric.  This  is  certainly  the  appearance  con¬ 
veyed  by  our  photographs  of  streamers,  and  it  is  what  one  expects  for  the  con¬ 
centrated  field  and  current  around  the  tip  of  a  narrow  filament.  To  obtain  a 
more  complete,  two-dimensional,  solution,  one  would  almost  certainly  use  a 
series  expansion  with  a  spherically  symmetric  current  as  the  leading  term. 

One  expects  this  to  be  dominant  in  the  region  so  close  to  the  streamer  that 
E  and  j  are  well  in  excess  of  the  background  field  and  current  of  the  unper¬ 
turbed  discharge.  Since  this  is  also  the  region  where  gas  heating,  which 
causes  streamer  growth,  is  concentrated,  keeping  only  the  radial  current 
should  give  a  reasonable  estimate  of  the  growth  rate. 

Assuming  that  the  current,  field,  and  gas  flow  are  radial  not  just  along 

e  =  0  but,  for  some  region  around  that  line,  one  sees  from  Faraday's  law 

Eq.  (69)  that  the  field  E  must  be  a  function  only  of  n,  and  from  charge 

n  2 

conservation  Eq.  (70)  that  the  current  density  must  vary  as  1/n  . 


Here  is  just  a  constant,  independent  of  both  time  and  position.  From 
Ohm's  law,  one  has 
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(73) 


Since  Euler's  equation  just  specifies  the  negligible  variations  in  the  pres¬ 
sure,  there  remain  only  two  equations:  the  continuity  equation  for  the  gas, 
which  now  has  the  simpler  form. 


-  I?  ^  ^  f;  <->'  0  =  “  <”> 

n 

and  the  enthalpy  Eq.  (68).  Since  j  is  known,  and  E  is  determined  by  a,  which 
is  determined  by  the  temperature  of  the  gas,  we  have  left  two  equations.  (74) 
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and  (68)  for  the  two  unknown  functions,  the  gas  velocity  v  (n)  and  the 

ff|0 

tenperature,  T  (n)  which,  at  constant  pressure,  p  ,  determines  the  density, 

0 

Internal  enthalpy,  and.  with  j,  the  conductivity  of  the  gas  at  points  around 
the  streamer. 

Equation  (74)  can  be  rewritten  as 

<''no  "  ■  ^*^0^]  *  -  YPn  On  +  Zc^) 

Integrating  along  e  ==  0  from  the  edge  of  the  core  n|^  to  a  point  n  In  the 
halo  or  beyond,  one  has 
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Here,  Is  the  density  of  the  gas  In  the  core  of  the  streamer,  at  points 
n  <  4nd  we  have  used  the  fact  that  In  the  core  heat  Is  radiated  or  con- 

“  U 

ducted  away  and  the  gas  flow  v  {n  -  n.)  goes  to  zero*  Remember  that.  In 

T|0  D 

general,  p  and  v  are  both  functions  of  the  radial  distance,  n-  Now  notice 
no 

that,  since  p^  Is  essentially  constant  within  the  arc-like  core,  because  the 
temperature  is  roughly  constant  there,  the  second  term  on  the  left  Is  Just  the 
Integral  from  zero  to  So  one  can  write  the  equation  as 
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The  assumption  of  spherical  symmetry  Is  not  expected  to  be  very  accurate 
within  the  core,  which  Is  a  cylindrical  structure,  although  the  cap  of  the 
core  may  have  a  somewhat  spherical  form.  But  formally  one  can  still  carry  the 
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Integral  on  1n  to  q  >  0  to  obtain  this  simpler  expression,  which  should  only 
be  used  for  points  q  outside  q^. 

Combining  this  result  with  Eq.  (68).  the  equation  for  the  enthalpy,  or 
gas  heating,  gives 


n 


E,  j. 


Here  the  gradient  of  H,  the  enthalpy,  has  been  written  as  a  total  derivative, 
since  everything  varies  only  with  radius  In  the  spherically  symmetric  solu¬ 
tion.  This  derivative  can  be  written  as 


dH  _  dH  iL 
dq  dT  dq 


because,  at  constant  pressure,  H  depends  only  upon  T,  provided  that  the 
kinetic  energy  of  gas  motion  Is  negligible,  which  It  1s,  to  good  approximation 
for  this  subsonic  flow.  Then  the  equation  becomes 
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since  p,  H,  and  a  are  all  known  functions  of  the  temperature,  T  (o  also  In¬ 
volves  the  current,  but  that  Is  known,  too)  and  t>  J  t  c  ,  and  q  (which 

0  0  c 

Is  needed  to  start  a  solution)  are  Just  numbers,  we  have  now  reduced  the  prob¬ 
lem  to  a  single  ordinary  differential  equation  for  a  single  unknown  function 
of  one  variable,  T(q).  This  equation  can  certainly  be  solved. 


s 
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5.0  DISCUSSION  OF  RESULTS 


Equation  (77)  can  be  solved  numerically  by  stepping  out  a  solution, 

T(n).  One  must  start  such  a  calculation  at  the  Inner  boundary,  since 
the  Integral  In  the  denominator  Involves  p  and.  hence,  T,  at  radii  smaller 
than  that  of  the  point  where  one  Is  calculating  dT/dn-  (Physically  this 
occurs  because  the  outward  radial  flow  of  gas  Is  caused  by  heating  Interior  to 
the  region  considered.)  So  one  starts  at  the  edge  of  the  core,  n^,  where 
the  temperature  Is  the  Ionization  temperature  of  the  gas,  and  steps  out¬ 
ward  In  n.  Since  y  Is  not  known,  this  must  be  done  for  a  whole  range  of  y 
values.  The  choice  of  y.  which  causes  T  to  approach  the  ambient  temperature 
of  the  background  discharge  at  large  distances  from  the  streamer.  Is  the  pre¬ 
dicted  growth  rate,  the  solution  to  the  problem. 

To  do  this  calculation  one  needs  j  ,  the  constant  that  specifies  the 

O 

streamer  current.  This  cannot  be  calculated  exactly  without  solving  the  whole 
two  dimensional  problem,  but  It  can  be  estimated.  As  discussed  earlier,  we 
expect  the  streamer  to  collect  approximately  that  current  that  would  otherwise 
have  been  drawn  by  the  area  of  the  cathode  that  Is  shadowed  by  the  streamer. 

I  *  A  V  C^  o  E  «  A  V  c^  o  E  e^^  (78) 

0  0  0  0  0 

Here,  as  before,  o  and  E.  are  the  conductivity  and  electric  field  of  the 
0  0 

unperturbed  discharge  and  A  Is  a  constant  near  unity.  In  the  absence  of  a 
two-dimensional  analysis,  A  should  be  taken  as  one,  unless  comparison  with 
data  shows  that  a  slightly  different  value  Is  more  accurate. 

Assuming  that  this  current  Is  collected  by  the  top  half  of  the  streamer 
halo,  one  has 


Combining  these  equations  gives 


Using  this  In  Eq.  (77). 
variable 


(81) 


Then  Eq.  (77)  becomes 


dx 


.2  2  p2 

_ 0  0 


4y  «  X^  ^  0  (3X'^  +  2X')  dX' 


(82) 


In  calculating  solutions,  this  equation,  not  Eq.  (77).  Is  the  most  convenient 
form  to  use. 

In  stepping  out  a  solution,  one  starts  at  the  Inner  boundary 

X  -  \  (83a) 

T  =  Tj  (83b) 

Here,  X^  Is  the  ratio  of  the  radius  of  the  tip  of  core,  exp  (yt).  to 
the  length  of  the  streamer,  c^  exp  (rt).  In  the  photographs,  the  ratio 
appears  to  be  of  the  order  of  0.1  or  less. 

To  do  the  numerical  calculation,  one  needs  to  know  the  core  radius.  The 
core  begins  where  the  temperature  reaches  T^  the  Ionization  temperature 
given  by  the  Saha  equation  for  the  discharge  gas  and  pressure.  At  this  point, 
the  field  will  drop  below  the  breakdown  field,  and  the  Ionization  will  be 


The  ratio  of  core  radius  to  streamer  length  Is,  obviously  enough.  Just  the 
square  root  of  the  ratio  of  discharge  current  density  to  core  current  den¬ 
sity.  Since  the  seen  In  the  photographs  was  about  0.1  and  the  discharge 

current  density  was  an  ampere  per  square  centimeter  or  more,  the  filamentary 

2 

core  of  the  streamer  must  have  carried  a  current  of  100  A/cm  or  more.  The 
reason  the  current  Is  this  high  Is  straightforward.  The  arc  behind  a  streamer 
develops  so  fast  that  there  Is  not  time  for  thermal  diffusion  to  spread  the 
current  very  much. 

There  Is  one  more  equation  or,  rather,  boundary  condition,  that  must 
also  be  considered.  This  Is  the  prescription  that,  since  the  streamer  core  Is 


a  good  conductor,  the  tip,  at  x^.  Is  essentially  at  cathode  potential. 

Hence  the  concentrated  electric  field  around  the  tip  must  drop  a  voltage  equal 

to  E  c,  the  voltage  that  would  have  been  dropped  over  the  streamer  length 
0 

In  the  background  discharge.  So  we  have. 


This  reduces  to 


2  dx 

ok  "J  .2 
0  \  X  a 

c 


(85) 


Notice  that,  since  this  Integral  contains  o.  the  spatially  varying  conduc¬ 
tivity.  which  depends  upon  the  gas  density,  which  depends  upon  T.  this  equa¬ 
tion  can  only  be  evaluated  after  one  has  solved  Eq.(82)  for  the  temperature 
profile,  T(X). 

One  cannot  be  completely  sure  until  one  has  generated  solutions  and 
studied  their  form,  and  the  present  program  did  not  get  that  far,  but  It 
appears  that  Eqs.(82)  and  (85)  together  are  enough  to  allow  calculation  of 
both  T  and  x^.  That  Is,  for  each  choice  of  y  and  X^,  one  obtains,  from 
the  resulting  temperature  profile,  T(X),  two  numbers,  the  temperature  at  x  =  1 
and  the  voltage  dropped  around  the  streamer,  which,  in  these  terms.  Is  the 
Integral  In  Eq.(85).  The  solution  to  the  problem  Is  that  pair,  y  and  x^, 
which  makes  T(1)  equal  to  the  ambient  temperature  and  satisfies  Eq.(85). 

There  Is  a  way  to  look  at  the  physics  of  this,  which  may  make  the  math¬ 
ematics  more  understandable.  The  streamer  grows  by  heating  the  gas  to  the 
conduction  temperature,  T^.  It  does  this  by  concentrating  the  current, 
which,  of  course,  gives  a  much  higher  heating  ahead  of  the  streamer  than  else¬ 
where  In  the  discharge.  One  might  then  wonder  why  the  streamer  does  not  *neck 
down*  even  more,  concentrating  the  current  Into  an  extremely  tiny  filament. 
Then  It  would  have  to  heat  even  less  gas.  and  the  current  density. 
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and  hence  the  rate  of  heating  at  the  tip  would  be  even  higher,  so  the  streamer 
could  grow  even  faster.  What  Eq.(85)  says  Is  that  this  does  not  happen  because 
the  streamer  must  concentrate  the  current,  using  only  the  available  voltage. 

If  the  streamer  made  a  narrower  core  and  grew  faster.  It  would  be  growing 
through  colder  gas,  which  would  have  a  lower  conductivity.  To  concentrate  the 
current  would  then  require  a  stronger  electric  field,  which  cannot  be  sup¬ 
ported  with  the  available  voltage.  The  point  Is  that  the  streamer  must  not 
only  heat  the  gas  that  eventually  becomes  part  of  the  core  but  It  must  also 
heat  all  the  gas  around  the  tip  enough  to  Increase  the  conductivity  of  the 
whole  surrounding  region,  so  that  the  current  can  all  be  funneled  Into  the  tip 
of  the  core  without  requiring  a  larger  voltage  drop  than  Is  available.  By 
considering  both  requirements.  It  appears  possible  to  calculate  both  the 
streamer  growth  rate  and  the  radius  of  the  tip  of  the  core. 

The  remaining  Inputs  needed,  namely  the  functions  p(T)  and  H(T)  and  the 
electrical  conductivity  a,  must  be  known  for  the  discharge  gas  mixture.  The 
density  and  enthalpy  are  standard  thermodynamic  functions.  The  conductivity, 
as  explained  earlier.  Is  given  In  the  halo  by  the  assumption  that  the  field  Is 
at  the  breakdown  threshold,  and  beyond  the  halo  by  the  e-beam-produced  elec¬ 
tron  density,  n  ,  and  the  mobility,  u  »  which,  of  course,  depends  upon  the 

*  B 

gas  density.  Since  one  does  not  a  priori  know  which  form  to  use,  the  proce¬ 
dure  Is  to  calculate  both,  then  use  whichever  conductivity  Is  higher. 

Finally,  It  Is  useful  to  consider  the  pressure-dependence  of  the 
result.  The  conclusion  of  the  scaling  analysis  In  Section  Z.O  was  that  one 
wants  to  go  to  higher  discharge  pressure.  Assume  that  the  pressure  Is  In¬ 
creased  by  a  factor  k.  At  constant  E/n,  the  field  E^  will  Increase  by  the 
same  factor.  If  the  e-beam  current  Is  Increased  proportionally,  the  electron 
density  In  a  recombination  dominated  discharge  will  also  be  Increased  by  a 

factor  k.  This  leaves  the  electrical  conductivity  o  unchanged,  since  the 

0 

mobility  Is  reduced  by  k.  So  the  discharge  current  Is  Increased  by  k,  keeping 
the  current  gain  constant. 

In  the  halo  of  a  streamer,  the  current  would  also  be  higher  by  a  factor 
of  k,  but  so  Is  the  breakdown  field.  Hence  the  conductivity  o  Is  unchanged 
(If  the  temperature  Is  similar).  From  Eq.  (82),  one  sees  that,  to  keep  the 
expression  similar,  one  would  have  to  Increase  y  by  a  factor  of  k  also.  The 
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unknown  here  Is  the  change  In  the  streamer  core  diameter.  But  even  al¬ 
lowing  for  this,  the  effect  on  the  equation  Is  worrysome. 

Consider  now  a  second  case,  where  the  e-beam  current  Is  left  at  the  same 
value  as  In  the  lower  pressure  discharge.  Then  the  electron  density  only 
Increases  as  the  conductivity,  0^,  decreases  by  and  the  discharge 

current  density,  E^,  and  hence  the  current  gains  of  the  system  are  In¬ 
creased  by  a  factor  yk.  Here  the  halo  conductivity  Is  reduced  by 
because  the  breakdown  field  Is  still  Increased  by  k,  but  the  current  only 
by  ^/ir for  streamers  of  similar  geometry.  Now,  to  keep  the  equation  the 
same,  one  must  again  Increase  y,  but  only  by  ^X.  The  unknown  factor  Is  the 
same:  the  change.  If  any.  In  X^. 

It  appears  that  this  higher  pressure  discharge  would  probably  arc 
sooner.  Basically,  what  the  equation  shows  Is  that.  If  the  discharge  specific 
power  Is  Increased,  the  streamer  growth  rate  should  be  higher.  That  was  true 
In  each  of  the  above  examples  since  the  E  field  was  Increased  by  the  same 
factor  as  the  pressure,  and  the  current  density  was  also  raised,  so  the  dissi¬ 
pation,  In  watts  per  gram  of  gas,  was  higher. 

This  Is  not  a  firm  conclusion,  since  the  change  In  x^  Is  not  known, 

and  we  only  considered  the  effect  In  the  halo,  not  the  lower  field  region 

beyond.  Indeed,  as  discussed  earlier,  there  Is  evidence  that  a  switch  type 

discharge  has  an  arcing  limit  that  depends  more  upon  field  than  power.  Still, 

2 

It  Is  plausible  that  If  E/p  Is  kept  constant  while  Is  raised, 

the  allowable  pulse  length  will  be  reduced. 

The  remaining  possibility  Is  to  reduce  E/p  somewhat  as  the  pressure  Is 
raised.  That  Is,  one  can  Increase  the  voltage  by  a  smaller  factor  than  k. 

Then  one  could  Increase  the  discharge  current  density  without  a  proportional 
Increase  In  specific  power.  In  terms  of  the  foregoing  analysis,  the  effect  of 
a  decrease  In  E/p  Is  to  decrease  the  size  of  the  high  field  halo  around  the 
tip  of  a  streamer.  The  calculations  that  we  have  made  strongly  suggest  that 
one  will  want  to  do  this. 

The  scaling  analysis  shows  that  this  technology  should  go  In  the  direc¬ 
tion  of  higher  pressure  discharges.  The  streamer  analysis  shows  that  It  will 
probably  not  be  feasible  to  Increase  the  voltage  by  as  large  a  factor  as  the 
pressure. 
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